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An  experimental  study  of  the  interaction  of  an  underwater  turbulent  round  jet  with  the  free 
surface  was  conducted.  Flow  visualization,  surface  curvature  measurements  and  hot  film 
velocity  measurements  were  used  to  study  this  flow.  It  is  shown  that  surface  waves  are 
generated  by  the  large  scale  vortical  structures  in  the  jet  flow  as  they  approach  the  free 
surface.  These  waves  propagate  at  an  angle  with  respect  to  the  flow  direction.  The 
propagation  angle  increases  as  the  strength  of  the  interaction  is  increased  by  increasing  the 
momentum  flux  of  the  jet  or  reducing  the  distance  of  the  jet  to  the  free  surface  or  both. 
Propagation  of  these  waves  in  the  flow  direction  is  suppressed  by  the  surface  current 
produced  by  the  jet.  Far  downstream  Lie  surface  motions  are  caused  by  the  large  scale 
vortical  structures  interacting  directly  with  the  surface.  The  fundamental  scaling  parameters 
of  the  firee-surface  jet  have  been  determined.  The  velocity  scale  is  the  velocity  obtained 
from  the  combination  of  jet  momentum,  density  and  depth  of  the  jet  and  the  length  scale  is 
the  distance  of  the  jet  to  the  free  surface.  It  is  shown  that  the  centerline  velocity  decay 
when  scaled  with  these  parameters  collapses  to  a  universal  curve  for  different  depths  of  the 
jet  The  asymptotic  decay  in  the  far  field  is  reduced  by  a  factor  of  2^^  compared  to  the  free 
jet  due  to  the  continement  by  the  free  surface.  The  growth  rate  of  the  free-surface  jet  is 
found  in  good  agreement  with  the  free  jet  However  the  eccentricity  of  the  jet  cross  section 
caused  by  the  displacement  of  the  jet  centerline  persists  for  large  distances  downstream, 
beyond  .40  times  the  initial  depth  of  the  jet  centerline.  Measurements  are  also  reported  on 
the  flow  field  of  a  jet  moving  parallel  to  a  solid  surface.  These  results  are  compared  with 
the  results  for  the  free-surface  jet.  It  is  found  that  the  solid  wall  jet  flow  is  fundamentally 
different  from  the  free-surface  jet  flow.  The  growth  rate  of  the  wall  jet  in  the  direction 
parallel  to  the  surface  is  approximately  3.9  times  the  growth  rate  of  the  free-surface  jet. 
The  growth  rate  of  the  wall  jet  in  the  direction  perpendicular  to  the  surface  is  half  the  free- 
surface  jet  growth  rate.  This  is  the  result  of  the  different  dynamics  of  vorticity  on  the  f%e 
surface  compared  to  a  solid  wall.  The  skin  friction  at  the  solid  wall  and  increased  growth 
rate  combined  to  give  a  different  maximum  velocity  decay  rate  compared  to  the  free-surface 
jet 
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CHAPTER  I 
INTRODUCTION 

The  results  of  an  experimental  investigation  of  the  interaction  of  a  submerged 
turbulent  jet  moving  parallel  to  the  free  surface  with  the  free  surface  are  presented.  The 
flow  geometry  and  parameters  are  shown  in  Figure  1.1.  One  motivation  for  this  study  is  to 
obtain  a  better  fundamental  understanding  of  the  nature  of  the  free  surface  waves  and 
motions  caused  by  the  interaction  of  turbulent  shear  flows  with  the  surface.  These  flow 
processes  are  an  important  part  of  the  free  surface  signature  of  the  turbulent  wake  behind  a 
ship.  Aerial  and  space  photographs  of  the  sea  surface  show  distinct  features  which  persist 
for  very  large  distances  behind  the  ship  in  the  viscous  wake  region  (Munk  et  al.  1987). 

The  turbulent  jet/free-surface  interaction  is  one  of  the  simplest  flow  configurations  which 
incorporates  many  of  the  vortical  interactions  encountered  in  the  turbulent  ship  wake 
problem. 

An  early  experimental  investigation  of  the  interaction  of  a  submerged  jet  with  the 
free  surface  was  conducted  by  Evans  (1955).  He  demonstrated  the  calming  effect  on 
surface  waves  caused  by  the  surface  currents  produced  by  the  jet.  Evans  study  did  not 
consider  the  details  of  the  turbulent  flow  structure.  He  showed  that  when  the  waves  and 
currents  move  in  the  same  direction  the  wave  amplitude  is  decreased.  A  theoretical  analysis 
of  this  phenomenon  by  Taylor  (1955)  provides  an  explanation  of  these  results  and  further 
show  that  in  this  case  the  wavelength  of  the  surface  waves  is  increased.  The  effect  of 
ntmuniform  steady  surface  currents  was  also  investigated  by  Longuet-Higgins  and  Stewart 
(1961)  theoretically  (See  also  Phillips, 1966).  They  consider  different  types  of  nonuniform 
surface  currents.  They  show  that  if  the  current  is  in  the  same  direction  as  the  direction  of 
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propagation  of  the  waves  the  amplitude  of  the  waves  decreases.  On  the  other  hand  if  the 
waves  propagates  into  a  region  where  the  surface  current  is  in  the  opposite  direction  as  the 
propagation  direction,  the  amplitude  of  the  waves  increases. 

Rajaratnam  and  Humphries  (1984)  studied  the  mean  flow  characteristics  of  free 
surface  jets  when  the  free  surface  is  located  at  the  edge  of  the  jet  nozzle,  i.e.  h/d=0.5.  In 
their  investigation  they  did  not  study  the  free  surface  motion  caused  by  the  jei/free-surface 
interaction.  However  they  reported  a  reduction  of  the  mean  velocity  near  the  surface  at 
high  Froude  numbers  which  was  attributed  to  surface  wave  generation.  For  circular 
surface  jets  they  confirmed  the  same  scaling  for  the  maximum  mean  velocity  decay  as  for 
the  free  jet  or  wall  jet  measurements  of  Rajaramam  and  Pani  (1974).  Self-similarity  was 
found  for  the  mean  velocity  profiles.  The  growth  rate  in  the  direction  perpendicular  to  the 
free  surface  was  found  equal  to  the  wall  jet  growth  rate  while  the  growth  rate  in  the 
direction  parallel  to  the  free  surface  was  found  to  be  approximately  half  of  the  wall  jet 
growth  rate.  Rajaratnam  and  Humphries  (1984)  and  more  recently  Ramberg  et  al.  (1989) 
have  studied  two-dimensional  frce-surface  jets.  They  also  pointed  out  the  similarity  with  a 
wall  jet  in  this  case.  Ramberg  et  al,  noted  the  pervasive  effects  of  jet  confinement  in  their 
tank.  These  confinement  effects  have  been  studied  by  Kotsovinos  (1976, 1978).  These 
effects  can  lead  to  break  down  of  the  similarity  scaling  laws  because  of  the  momentum  flux 
associated  with  the  entrained  fluid.  Novikov  (1988)  discussed  the  mean  surface 
deformation  of  the  free  surface  caused  by  an  underwater  jet  He  did  not  consider  the  effect 
of  the  velocity  fluctuations. 

The  significance  of  the  compariscm  between  the  wall  jet  and  the  free-surface  jet  is 
not  immediately  obvious.  The  main  similarity  between  these  flows  is  the  confinement  by 
the  surface  which  limits  the  flow  at  the  free  surface  or  at  the  wall.  However,  in  a  wall  jet 
the  no  slip  boundary  condition  requires  the  velocity  to  be  zero  at  the  wall  which  results  in  a 
boundary-layer-like  behaviour  near  the  solid  surface.  On  the  other  hand  at  a  free  surface 
the  velocity  can  be  different  from  zero.  The  viscous  boundary  condition  at  a  free  surface 


requires  continuity  of  normal  and  shear  stresses  at  the  water  surface  (Batchelor,  1967 
Sec.3.3).  For  a  flat  free  surface  and  neglecting  the  shear  stress  of  the  air  this  condition 
implies  zero  stress  at  the  surface.  Although  a  free  surface  boundary  layer  will  form  to 
satisfy  the  nonsteady  viscous  boundary  condition  at  the  free  surface  and  possible  capillary 
effects,  the  dynamics  of  this  free-surface  boundry  layer  is  fundamentally  different  than  the 
boundary  layer  at  the  solid  surface. 

Davis  and  Winarto  (1980)  studied  the  interaction  of  a  jet  with  a  solid  surface  for 
different  distances  between  the  jet  and  the  solid  surface.  This  geometry  is  the  same  as  the 
one  studied  in  the  present  investigation  except  for  the  fact  that  the  free  surface  is  replaced 
by  a  solid  surface.  The  results  of  Davis  and  Winarto  are  in  disagreement  with  the  results  of 
Rajaratnam  and  Pani  (1974).  Davis  and  Winarto  (1980)  measured  a  somewhat  higher 
growth  rate  in  the  direction  perpendicular  to  the  solid  wall  and  somewhat  lower  growth  rate 
in  the  direction  parallel  to  the  solid  wall  compared  to  the  the  results  of  Rajaratnam  and  Pani. 
This  work  as  well  as  other  investigations  of  the  solid  wall  jet  for  different  jet  exit 
geometries  by  Sforza  and  Herbst  (1970),  Newman  et  al.  (1972),  Chandrasekhara  and 
Bandyopadhyay  (1975)  was  reviewed  by  Launder  and  Rodi  (1981, 1983).  Launder  and 
Rodi  (1981)  noted  that  the  high  rate  of  spanwise  spreading  of  the  wall  jet  may  be  caused  by 
a  secondary  motion  in  the  cross-sectional  plane.  This  is  attributed  to  the  generation  of 
streamwise  vorticity  by  bending  of  vortex  lines  near  the  wall. 

In  contrast  to  the  turbulent-jet/free-surface  interaction  problem  the  flow 
characteristics  of  a  turbulent  free  jet,  i.e.  in  the  absence  of  the  free  surface,  has  been  the 
subject  of  iTUUiy  investigations.  It  seems  apnrqrriate  to  discuss  the  underwater  flow 
characteristics  in  the  turbulent  jet/free-surface  interaction  in  terms  of  the  flow  characteristics 
of  the  free  jet  problem.  The  flow  field  can  be  divided  into  three  regions:  the  near  field,  the 
transition  region  and  the  far  field  (Abramovitz  1963).  The  mean  and  turbulent  flow 
properties  in  the  far  field  of  a  turbulent  jet  have  been  measured  by  several  investigators  and 
nx)st  extensively  by  Wygnanski  and  Fiedler  (1969).  In  the  far  field  the  maximum  mean 
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velocity  at  each  cross-section  decays  like  1/x,  where  x  is  the  distance  to  the  jet  exit  plane, 
and  the  width  of  the  velocity  profile  is  proportional  to  x.  This  scaling  implies  that  the  flux 
of  downstream  momentum  is  constant  along  the  axis  of  the  flow.  As  mentioned  earlier  for 
the  two-dimensional  jet,  Kotsovinos  (1978)  and  Schneider  (1985)  have  investigated  the 
effect  of  the  momentum  flux  associated  with  the  entrained  fluid.  They  show  a  continuous 
reduction  of  the  momentum  flux  with  downstream  distance  due  to  this  effect 

The  turbulent  structure  of  the  free  jet  has  been  the  subject  of  many  investigations. 
In  the  near  field  the  potential  core  region  is  surrounded  by  a  turbulent  shear  flow  region. 
For  an  initially  laminar  shear  layer,  instability  waves  grow  exponentially  and  result  in  the 
formation  of  vortex-ring-like  structures  which  interact  with  each  other  by  amalgamation  as 
they  grow  downstream.  At  the  end  of  the  potential  core  in  the  transition  region  vorticity 
within  these  structures  reaches  the  centerline  and  associated  with  this  process  there  is  a 
characteristic  frequency  or  preferred  mode  of  the  jet  (Browand  and  Laufer,  1975,  Yule, 
1978).  In  the  far  field  of  the  jet  the  evidence  for  the  existence  of  large  scale  turbulent 
structures  is  convincing  (Tso,  Kovasznay  and  Hussain,1981  and  Dimotakis  etal.,  1983). 
There  are  however  a  number  of  unanswered  questions  regarding  the  topology  and 
dynamics  of  these  structures  in  the  far  field.  It  follows  that  an  important  aspect  of  the 
interaction  of  the  free-surface  with  the  turbulent  jet  is  to  determine  the  effect  of  the  free- 
surface  on  these  turbulent  structures  as  well  as  to  determine  the  role  of  these  structures  in 
the  dynamics  of  the  free-surface. 

In  the  present  investigation,  we  consider  the  interaction  of  a  circular  jet  with  the 
free  surface  fen*  various  depths  of  the  jet  below  the  surface.  The  main  objectives  of  this 
investigation  are: 

•  To  determine  the  surface  waves  and  motions  produced  by  the  interaction  of  the 
turbulent  jet  with  the  free  surface. 

•  To  determine  the  scaling  characteristics  of  the  underwater  turbulent  jet. 
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•  To  investigate  the  differences  and  similarides  between  the  free-surface  jet  and  the 
wall  jet  to  clarify  the  role  of  the  boundary  conditions  at  the  surface  in  these 
flows. 

The  presentation  of  the  results  is  organized  as  follows.  In  Chapter  n  the  flow 
facility  and  instrumentation  is  described.  In  Chapter  III  the  results  of  a  flow  visualization 
study,  surface  curvature  measurements  and  hot  film  velocity  measurements  for  the  free- 
surface  jet  and  the  wall  jet  are  presented.  In  Chapter  IV  these  results  are  discussed  and 
finally,  in  Chapter  V,  the  main  conclusions  of  the  investigation  are  summarized. 


CHAPTER  n 


EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 

n.l  Flow  Facility 

The  experiments  were  conducted  in  a  water  tank  facility  specifically  designed  for 
flow  visualization  and  Laser  diagnostics.  Some  of  the  criteria  considered  in  designing  this 
facility  are  the  ability  to  maintain  a  constant  water  temperature  and  quality,  as  well  as  the 
ability  to  maintain  a  constant  water  level  in  the  free  surface  tank.  Photographs  of  the 
facility  are  shown  in  Figure  2.1  and  Figure  2.2.  A  schematic  diagram  of  the  facility  is 
shown  in  Figure  2.3.  It  consists  of  a  free  surface  tank,  a  reservoir  tank,  a  jet  tank  and  the 
associated  piping  and  control  valves.  Tlic  free  surface  tank  was  made  of  glass  76.2  cm 
wide,  76.2  cm  high  and  167.6  cm  long.  The  frame  for  this  tank  was  constructed  from  steel 
angles  5.08  cm  x  5.08  cm  x  0.64  cm  and  steel  plates.  The  frame  was  bolted  to  a  double 
Unistrut'  structure  which  raised  it  80  cm  above  the  floor  for  optical  access  through  the 
bottom  surface.  The  steel  frame  was  isolated  from  the  glass  walls  by  0.32  cm  thick  open 
cell  sponge  rubber.  The  side  walls  were  made  of  1.27  cm  thick  glass.  The  bottom  surface 
was  made  in  two  sections.  A  glass  section  1.90  cm  x  147.3  cm  x  76.2  cm  and  a  PVC 
section  1.90  cm  X  20.3  cm  x  76.2  cm.  Drainage  and  water  intake  are  provided  by  a  5.08 
cm  diameter  hole  and  a  1.90  cm  hole,  respectively,  located  on  the  PVC  section  of  the 
bottom  face. 

Water  used  in  the  facility  was  stened  in  a  1800  liter  reservoir  tank.  The  external 
dimensions  of  the  reservoir  tank  are  122  cm  in  width,  61  cm  in  height  and  244  cm  in 
length.  The  tank  was  constructed  from  0.64  cm  thick  sheets  of  PVC.  To  further 
strengthen  the  tank,  it  was  braced  by  two  Unistrut'  structures  one  located  at  the  top  and  the 
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other  at  the  mid-height  of  the  tank.  The  tank  was  covered  with  0.64  cm  thick  PVC  sheet  to 
prevent  dust  particles  and  other  debris  from  contaminating  the  water. 

The  jet  flow  was  generated  by  a  jet  tank  located  inside  the  free  surface  tank.  A 
diagram  of  the  jet  tank  is  shown  in  Figure  2.4.  Water  entered  the  tank  from  above.  Large 
scale  flow  nonuniformities  associated  with  the  inlet  flow  were  removed  by  a  layer  of  foam 
5.08  cm  thick  located  between  two  Acrylic  sheets  with  square  meshes  of  1.27  cm  x  1.27 
cm  as  shown  in  Figure  2.4.  The  jet  exit  orifice  is  located  on  the  side  wall  at  the  mid-point 
between  the  foam  and  the  bottom  of  the  tank.  A  circular-arc-shaped  nozzle  with  a  radius 
equal  to  the  wall  thickness  provides  a  smooth  transition  from  the  side  wall  of  the  tank  to  the 
jet  exit  plane.  In  continuous  operation  the  large  area  ratio  of  the  contraction  helps  reduce 
the  turbulent  intensity  at  the  jet  exit  to  a  level  lower  than  0.5%.  Three  jet  tanks  were  used 
in  the  experiments  with  jet  exit  diameters  of  2.54  cm,  1.27  cm  and  0.64  cm  respectively. 
The  geometrical  dimensions  of  these  tanks  are  presented  in  Table  2. 1 .  The  jet  tank  for  the 
1.27  cm  exit  diameter  jet  was  construaed  from  1.27  cm  thick  Acrylic  sheets  and  was  used 
for  flow  visualization  and  surface  curvature  measurements.  The  other  two  were 
constructed  from  PVC  and  Acrylic  plate. 

The  velocity  profile  at  the  exit  plane  of  the  jet  was  uniform  except  for  a  boundary 
layer  region  at  the  surface.  The  momenmm  thickness  of  the  boundary  layer,  0,  was 
estimated  using  Thwaites  method.  The  result  for  the  three  jet  tanks  are  presented  in  Table 
2.2.  The  calculated  values  of  the  boundary  layer  Reynolds  number  based  on  the 
displacement  thickness,  5i,  suggests  that  the  boundary  layer  is  laminar  except  for  the 
higher  values  of  the  velocity  when  it  approaches  the  instability  point  ( See  for  example 
Schlichting,  1955.) 

Experiments  were  conducted  at  several  jet  exit  velocities,  Uc,  and  jet  exit  depths,  h. 
The  jet  exit  velocity  was  adjusted  with  suitable  needle  valves  located  in  the  facility's  control 
panel.  The  jet  exit  depth  was  kept  at  the  desired  value  by  maintaining  a  constant  level  in  the 
free-surface  tank  during  the  experiment  A  schematic  diagram  of  the  control  panel  and  the 
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associated  piping  is  shown  in  Figure  2.5.  In  addition  to  the  flow-control  needle  valves  the 
control  panel  includes  the  necessary  plumbing  for  various  house  keeping  chores. 

The  facility  was  operated  in  two  different  modes,  an  open  loop  configuration  and  a 
closed  loop  configuration.  The  open  loop  configuration  was  used  in  several  flow 
visualization  experiments  and  in  the  surface  curvature  measurements.  In  this  case  tap  water 
is  stored  in  the  reservoir  tank  and  pumped  through  the  facility  into  the  drain.  The  water 
level  in  the  free  surface  tank  was  kept  constant  by  means  of  a  stand-up  pipe  technique 
shown  schematically  in  Figure  2.6.  In  the  open  loop  operation,  the  PVC  discharge  valve  is 
closed  so  that  excess  water  above  the  desired  level  overflows  through  the  stand-up  pipe 
into  the  level  control  box  and  into  the  drain.  In  the  closed  loop  configuration  water  in  the 
firee-surface  tank  is  recirculated  through  the  jet  tank  using  the  pump.  The  closed  loop 
operation  was  used  in  velocity  measurements  to  maintain  a  constant  temperature  and  water 
quality  throughout  the  system  during  the  experiment.  The  maximum  flow  rate  that  can  be 
obtained  using  the  pump  in  this  facility  was  15  lit/min.  However,  operating  the  facility  at 
such  a  high  flow  rate  caused  excessive  recirculation  in  the  free  surface  tank.  Therefore,  the 
facility  was  operated  at  considerably  lower  flow  rates.  The  entrainment  velocity  was 
calculated  using  the  results  of  Ricou  &  Spalding  (1961).  The  value  of  the  recirculation 
velocity  for  Uc=2()0  cm/s,  d=0.64  and  x/d=40  was  estimated  to  be  0.23  cm/s  which  is 
apim^ximately  0.7%  of  the  local  maximum  mean  velocity.  The  maximum  jet  exit  velocity 
used  was  220  cm/s  in  the  0.64  cm  diameter  jet  tank. 

Experiments  were  also  conducted  on  the  interaction  of  the  je^  with  a  flat  solid 
surface.  A  schematic  diagram  of  the  experimental  apparatus  used  is  shown  in  Figure  2.7. 
The  apparatus  and  instrumentation  used  were  the  same  as  for  the  jet  /  free-surface 
interaction  experiments,  with  the  addition  of  a  solid  wall.  The  solid  wall  was  made  of 
0.64  cm  thick  PVC  plate  61  cm  long  and  55  cm  wide.  In  order  to  stiffen  the  plate,  PVC 
angles  5.08  cm  x  5.08  cm  x  0.64  cm  were  glued  on  the  edges  and  one  on  the  center  of  the 
plate.  The  surface  of  the  plate  was  smoothened  by  polishing  it  with  steel  wool.  The  plate 
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was  suspended  from  the  edges  of  the  free  surface  tank  by  four  1.27  cm  diameter  stainless 
steel  threaded  rods.  It  was  then  set  against  the  jet  tank  at  the  desired  height  below  the  jet 
nozzle.  The  small  gap  between  the  plate  and  the  jet  tank  was  covered  by  a  water-proof 
transparent  tape.  The  plate  was  made  flat  to  within  0.3  mm  over  its  total  length  of  60. 1  cm. 
In  these  experiments  the  free  surface  was  located  at  least  28  jet  exit  diameters  from  the 
centerline. 

II.2  Flow  Visualization 

Flow  visualization  of  the  free  surface  deformation  and  of  the  jet  fluid  was  obtained 
using  the  shadowgraph  technique.  A  schematic  diagram  of  the  optical  system  used  is 
shown  in  Figure  2.8.  A  collimated  beam  of  light  30  cm  in  diameter  was  formed  from  the 
output  beam  of  a  Copper  Vapor  Laser  2.2  cm  in  diameter  by  means  of  a  focussing  lens  and 
a  spherical  mirror.  A  30  cm  x  30  cm  first  surface  mirror  was  utilized  to  direct  the 
collimated  beam  of  light  through  the  bottom  glass  window  perpendicular  to  the  undisturbed 
free  surface.  Refraction  at  the  water-air  interface  causes  variations  in  beam  intensity  in  the 
regions  where  the  free  surface  is  not  perpendicular  to  the  beam.  The  resulting  shadow 
image  was  viewed  on  a  screen  located  above  and  close  to  the  free  surface.  The  pictures 
presented  here  were  obtained  by  photographing  the  image  on  the  screen  with  a  35  mm 
camera.  The  jet  fluid  was  simultaneously  visualized  by  using  1-2  *C  warmer  water  in  the 
jet  tank.  The  corresponding  change  in  the  index  of  refraction  results  in  the  light  intensity 
variation  typical  of  shadowgraph  images.  Motion  pictures  of  this  flow  were  also  obtained 
using  a  Hycam  Model  K20S4E  high  speed  camera  at  a  nominal  frame  rate  of  100  pictures 
per  scond. 

Laser  Induced  Fluorescence  (LIF)  was  used  to  visualize  cross-sections  of  the 
underwater  flow.  A  schematic  diagram  of  the  apparatus  used  for  LIF  flow  visualization  is 
shown  in  Figure  2.9.  A  fluorescent  dye  (Rhodamine  6G)  with  a  concentration  of  2  ppm 
was  homogeneously  mixed  with  water  in  the  jet  tank.  The  dyed  water  was  then  driven 


through  the  jet  nozzle  into  clear  water  in  the  free  surface  tank.  Excitation  of  the 
fluorescence  was  induced  by  a  plane  sheet  of  light  from  a  Copper  Vapor  Laser  formed  by  a 
combination  of  cylindrical  lenses  with  focal  lengths  of  -60mm,  100mm  and  150mm 
corresponding  to  lenses  LI,  L2  and  L3,  respectively,  in  Figure  2,9.  This  combination  of 
lenses  allowed  control  of  the  thickness  and  width  of  the  laser  sheet.  The  sheet  thickness 
was  1  -2  mm  over  the  region  of  observation.  Rhodamine  6G  dye  fluoresces  efficiently 
when  excited  by  the  green  (511  nm)  line  of  the  Copper  Vapor  Laser.  The  scattered  light, 
yellowish  in  color,  from  dye-containing  fluid  in  the  plane  of  the  laser  sheet  was  then 
recorded  by  a  35  mm  camera  positioned  perpendicular  to  the  sheet 


n.3  Surface  Curvature  Measurements 

Surface  curvature  measurements  were  conducted  to  quantify  the  various  surface 
features  observed  in  the  flow  visualization  study.  A  schematic  diagram  of  the  apparatus 
used  for  these  measurements  is  shown  in  Figure  2.10.  The  technique  used  is  based  on  the 
same  principle  as  the  shadowgraph  flow  visualization  technique.  In  this  case  a  photodiode 
(SKAN-A-MATIC  Model  SI  18  1/4)  was  used  to  measure  the  temporal  fluctuation  of  light 
intensity  at  a  fixed  point  on  the  surface.  The  diameter  of  the  photocell  was  5  mm  and  it  had 
a  nominal  frequency  response  of  300  KHz.  A  1  mm  diameter  circular  aperture  was 
utilized  to  improve  the  spatial  resolution  of  the  measurements.  In  the  experiments  the 
photodiode  was  positioned  0.64  cm  above  the  water  surface.  A  collimated  beam  of  light 
was  formed  using  the  514.5  nm  line  of  an  argon-ion  laser  (Lexel  Model  95)  as  the  light 
source.  Following  the  notation  in  Figure  2.1 1  the  light  intensity  at  the  photodiode  aperture 
from  a  two-dimensional  deformation  of  the  surface  t)  is  given  by 

lu')  =  1 1^1 

0  dx 

where  x'  is  the  location  of  the  photodiode  and  x  is  the  location  where  the  ray  captured  by 
the  photodiode  intersects  the  surface.  In  the  limit  of  small  surface  deformations  ^ 
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small  surface  slope  g;f«  1  and  small  surface  curvature  s — ^«1,  it  can  be  shown  that 

dx"^ 


the  light  intensity  measured  by  the  photodiode  is  given  by 


I(x’)  = 


-  1) 


d_^' 

dx2  / 


where  s  is  the  distance  from  the  photodiode  to  the  water  surface  and  nw  is  the  index  of 
refraction  of  water.  Thus  in  this  limit  the  light  intensity  measured  by  the  photodiode  is 
proportional  to  the  surface  curvature.  A  detailed  analysis  of  this  relationship  between  the 
light  intensity  measured  by  the  photodiode  and  the  surface  deformation  is  presented  in 
Appendix  A. 

Figure  2.12  is  a  block  diagram  showing  the  data  acquisition  and  storage  system. 
The  signal  from  the  photocell  was  DC  shifted,  amplified  and  filtered  using  a  Tektronix  AM 
502  differential  amplifier.  The  analog  ouqjut  of  the  filter  was  then  digitized  using  a  Lccroy 
Model  8210  Transient  Digitizer.  Typical  sampling  rate  was  200  Hz.  The  digitized  output 
was  then  stored  on  permanent  files  using  an  IBM  CS9000  mini  computer. 


n.4  Hot  Film  Velocity  Measurements 

Velocity  data  were  obtained  using  a  constant  temperature  hot  film  anemometer.  A 
standard  TSI  quartz  coated  cylindrical  sensor  with  a  frequency  response  of  approximately 
80  KHz  was  used  in  the  measurements.  The  sensor  length  was  0.51  mm  and  the  diameter 
25  ^m.  The  sensor  axis  was  positioned  perpendicular  to  the  flow  direction  and  parallel  to 
the  free  surface.  The  hot  film  was  operated  at  the  oveiheat  ratio  of  1.09.  Figure  2. 13  is  a 
block  diagram  of  the  electronics  and  data  acquisition  system.  A  TSI  Model  1750  constant 
temperature  anemometer  box  in  conjunction  with  a  TSI  Model  1056  Decade  Resistance 
box,  with  a  resolution  of  0.01  was  used  to  obtain  the  data.  The  output  of  the 
anemometer  was  DC  shifted  and  amplified  using  a  Tektronix  AM501  operational  amplifier 
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wired  as  a  differential  amplifier  with  a  gain  of  2.6.  The  output  of  the  differential  amplifier 
was  digitized  using  a  Lecroy  Model  8210  Transient  Digitizer.  Typical  sampling  rates  used 
were  between  200-800  Hz.  The  digitized  output  was  then  stored  on  permanent  files  using 
an  IBM  CS9000  computer. 

A  schematic  diagram  of  the  apparatus  used  for  velocity  measurements  is  shown  in 
Figure  2.14.  The  hot  film  probe  was  mounted  on  a  computer  controlled  traverse 
mechanism.  The  traversing  gears  in  the  axial  (x-axis)  and  transverse  (y-axis)  directions 
were  actuated  by  two  Sigma  Model  18-1408  D40-F  stepping  motors  with  a  50:1  gear  ratio 
(Willmarth,  1977).  A  Transicoil  Model  U-2 17094  Gearhead  Stepper  Motor  with  a  gear 
ratio  of  60:1  was  used  to  move  the  probe  in  the  direction  perpendicular  to  the  free  surface. 
The  stepping  motor  was  controlled  by  the  IBM  CS9000  computer  using  a  parallel  port.  An 
interface  circuit  was  designed  and  built  to  conven  the  parallel  pon  output  to  the  required 
control  signals  for  the  stepping  motor.  It  provided  full  three-dimensional  positioning  of  the 
probe.  The  resolution  of  the  traverse  motion  along  the  axial  and  transverse  directions  was 
25  pm  and  in  the  vertical  direction  was  6.86  pm. 

The  hot  film  probe  was  calibrated  in  the  same  free  surface  tank  facility  used  in  the 
experiments  by  locating  it  at  the  jet  exit  and  varying  the  jet  exit  velocity.  The  jet  exit 
velocity  was  measured  by  a  manometer  which  was  constructed  using  a  precision  Vernier 
Caliper,  0.64  cm  diameter  glass  tubing  and  0.95  cm  diameter  Tygon  tubing.  A  schematic 
diagram  of  the  pressure  manometer  is  shown  in  Figure  2.15.  The  height  of  the  column  of 
water  cwresponding  to  the  total  pressure  Pt  and  the  static  pressure  at  the  exit  of  the  jet 
were  measured  using  the  manometer.  The  difference  between  the  total  and  static 
pressures,  measured  in  cm  of  water  is  proportional  to  the  square  of  the  jet  exit  velocity. 

For  each  jet  exit  velocity,  the  output  voltage  of  the  amplifier  and  the  exit  velocity  obtained 
from  the  manometer  were  recorded.  A  typical  calibration  curve  is  shown  in  Figure  2. 1 6. 
The  solid  line  tm  this  plot  is  a  fourth  order  polynomial  fit  through  the  measured  data 
obtained  using  the  least  squares  method.  The  triangles  in  this  figure  represent  the  measured 


13 


data  points.  The  error  between  the  measured  and  calculated  values  is  less  than  1%  of  the 
maximum  jet  exit  velocity  Ue,„.  The  calibradon  curve  was  incorporated  into  the  Data 

Acquisition  software  and  used  to  obtain  the  velocity  for  each  measured  value  of  the  hot  film 
output  voltage. 

In  order  to  minimize  probe  interference  with  the  free  surface  at  close  distances  to  the 
free  surface,  the  probe  was  tilted  backward  by  4  degrees  (See  Figure  2.14).  The  velocity 
calibration  was  performed  with  the  probe  in  the  tilted  position.  The  closest  point  to  the 
water/air  interface  for  which  the  measurements  were  performed  was  at  0.64  mm  from  the 
surface.  At  this  point  the  hot  film  sensor  was  intermittently  in  and  out  of  the  water,  which 
resulted  in  sharp  voltage  spikes  on  the  time  trace  of  the  velocity  signal.  The  velocity  data 
corresponding  to  the  points  close  to  the  surface  were  examined  and  those  contaminated 
with  sharp  spikes  were  discarded. 

For  the  wall  jet  velocity  measuren%nts,  the  probe  support  was  modified  in  order  to 
gain  access  to  the  jet  exit  orifice  and  to  conduct  measurements  as  close  as  possible  to  the 
wall.  Referring  to  Figure  2.7,  the  90  degrees  angle  adapter  was  removed  from  the  probe 
support.  The  probe  was  mounted  directly  on  the  probe  support  which  was  tilted  backward 
at  an  angle  of  9  degrees  as  shown  in  Figure  2.7.  The  closest  point  to  the  flat  plate  for 
which  data  could  be  obtained  was  at  0.25  mm  above  the  plate.  The  probe  was  calibrated  in 
this  configuration  by  position  it  at  the  exit  of  the  jet 

In  the  velocity  measurements  particular  emphasis  was  placed  on  insuring  a  constant 
temperature  and  quality  of  the  water  in  the  tanks  to  minimize  drift  of  the  velocity  signal. 

The  water  temperature  and  the  probe's  cold  resistance  were  monitored  during  the 
measurements.  The  water  temperature  was  maintained  constant  to  within  0.5  °C.  In  order 
to  improve  the  water  quality,  the  tap  water  was  passed  through  a  nominal  ft  pm  sediment 
rernoval  filter  and  a  water  softener  filter  before  filling  the  reservoir  tank.  To  inhibit  the 
growth  of  one-ceiled  animal  and  plant  life.  Chlorine  with  the  concentration  of  3  ppm  was 
added  to  the  water  in  the  reservoir  tank.  Hie  reservoir  tank  was  covered  with  a  sheet  of 
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gray  PVC  0.64  cm  thick  to  prevent  light,  dust  and  dirt  particles  from  falling  into  the  clean 
water.  The  water  stored  in  the  reservoir  tank  was  used  to  fill  the  fice  surface  tank  to  the 
desired  height.  The  jet  was  produced  by  recirculating  the  clean  treated  water  in  the  free 
surface  tank  using  the  pump.  To  insure  that  no  significant  drift  in  the  hot  film 
instrumentation  occurred,  the  probe  was  moved  to  the  jet  exit  orifice  at  the  end  of  each 
traverse  and  the  jet  exit  velocity  measured.  If  a  significant  change  in  the  measured  velocity 
at  the  jet  exit  were  detected  the  measurements  were  discharged. 


CHAPTER  in 


RESULTS 


m.l  Flow  Visualization 

Figures  3.1  and  3.2  show  LIF  photographs  of  the  jet  development  and  interacdon 
of  the  jet  with  the  free  surface.  In  both  figures  the  jet  flow  is  from  left  to  right  and  the  jet 
centerline  was  positioned  one  diameter  below  the  surface.  Figure  3.1  shows  a  cross- 
section  view  of  the  underwater  jet  through  the  symmetry  plane  at  a  Reynolds  number  of 
approximately  2.5x1  (P.  The  vortex  ring  like  structures  in  the  potential  core  of  the  jet 
(0<x/d<4)  are  clearly  observed  in  this  image.  As  they  cmivect  downstream  these  vortical 
structures  interact  with  each  other  and  with  the  surface.  From  this  image  the  interaction  of 
the  jet  with  the  surface  begins  at  approximately  x/ds4.  Downstream  of  this  point  smaller 
scale  features  are  observed.  Hgure  3.2  shows  a  similar  cross-section  obtained  at  a 
Reynolds  number  of  8.9xl(P.  In  this  case  the  vcntical  structures  are  formed  with  a 
wavelength  smaller  than  in  Hgure  3.1.  The  wavelength  of  the  vortical  structures  increases 
with  downstream  distance.  The  interaction  with  the  free  surface  starts  also  at 
approximately  x/d=4,  the  same  value  as  for  the  flow  in  Hgure  3.1.  The  wavelength  of  the 
stroctures  at  the  beginning  of  the  interaction  with  the  free  surface  is  approximately  the  same 
in  both  pictures. 

Typical  flow  visualization  shadowgraph  pictures  of  jet^ree-surface  interaction  are 
presented  in  Figures  3.3  to  3.8.  In  all  cases  the  flow  is  from  left  to  right.  The  jet  fluid  is 
visualized  by  a  slightly  lower  temperature  of  the  jet  tank  fluid.  Figures  3.3  to  3.5  were 
obtaiiwd  at  a  jet  exit  depth  corresponding  to  h/d=l.  In  these  cases  the  interaction  of  the  jet 
with  the  free  surface  occurs  in  the  near  Held  of  the  jet.  Figure  3.3  shows  the  interaction  at 
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low  jet  exit  velocity  as  indicated.  Near  the  jet  exit  the  vortex-ring-like  structures  can  be 
observed.  These  structures  interact  with  the  free  surface  and  produce  characteristic  surface 
deformations  which  travel  with  the  vortices.  The  light  intensity  pattern  along  the  jet 
centerline  in  Figure  3.3  indicates  a  depression  of  the  surface  on  top  of  the  core  of  the 
vortices  and  an  elevation  of  the  surface  between  the  vortices.  Further  downstream  the  jet 
vortical  fluid  interacts  with  the  surface.  In  this  region  the  shadowgraph  image  of  the  jet 
fluid  is  dominated  by  the  small  scale  features  in  the  flow.  A  remaiicable  feature  in  this 
photograph  are  the  two  dark  spots  located  at  x/da4.  The  dark  spots  are  associated  with 
vortex  lines  terminating  at  the  free  surface  (Berry  and  Hajnal  1983  and  Sterling  et  al 
1987).  They  are  generated  by  opening  of  vortex  lines  of  the  initially  submerged  coherent 
eddy.  This  type  of  interaction  has  been  investigated  by  Bernal  and  Kwon  (1989)  and 
Kwon  (1989)  for  axisymmetric  vonex  rings. 

The  flow  pictures  in  figures  3.4  and  3.5  were  obtained  at  jet  exit  velocities  of  35 
cm/s  and  50  cm/s  respectively.  These  pictures  show  typical  surface  wave  patterns 
produced  by  the  interaction  of  the  jet  with  the  free  surface.  The  flow  conditions  in  Figure 
3.4  are  slightly  above  the  values  for  which  surface  waves  are  first  observed.  The  vortex¬ 
ring-like  structures  moving  underneath  the  surface  start  deforming  the  surface  at  a  distance 
of  x/d»l,  and  it  is  not  before  a  distance  of  x/d=»4  for  which  waves  are  produced.  At  these 
flow  conditions  the  waves  are  produced  in  the  region  where  the  coherent  eddy  opens  up  in 
Figure  3.3.  The  resulting  wave  fronts  propagate  at  an  angle  of  39  degrees  with  respect  to 
the  downstream  direction.  Downstream  of  the  wave  generation  region  the  turbulence  in  the 
jet  interacts  with  the  free  surface.  Dark  spots  indicating  vortex  lines  terminating  at  the  free 
surface  are  observed  in  this  region.  At  a  higher  jet  exit  velocity  in  Figure  3.5  the  waves  are 
generated  over  a  larger  region  in  the  downstream  direction.  In  this  case  the  waves  are 
initially  formed  at  x/d«4  and  propagate  at  an  angle  of  60  degrees  relative  to  the  downstream 
direction.  However  the  wavefronts  curved  farther  downstream  suggesting  a  steepening  of 
the  direction  of  propagation.  Along  the  centerline  the  wavefronts  are  quite  well-defined  and 


are  nonnal  to  the  flow  direction.  This  picture  and  the  motion  pictures  of  the  flow  obtained 
at  the  same  flow  conditions  show  an  increase  of  the  wavelength  of  the  waves  as  they 
propagate  along  the  centerline. 

The  same  basic  features  observed  in  die  near  field  interaction  are  also  found  when 
the  interaction  occurs  in  the  far  field  of  the  jet  Figure  3.6  is  an  instantaneous  (l-2|is) 
spark  shadowgraph  picture  of  the  flow  field  at  a  Reynolds  number  of  1.27x10^  and 
h/d=3.5.  The  corresponding  jet  exit  diameter  was  d=1.27  cm.  The  shadowgraph  image  of 
the  flow  shows  small  scale  turbulent  structures  in  the  submerged  jet  fluid.  Although  more 
remarkable  is  the  deformation  of  the  surface  at  x/d»16  which  is  caused  by  the  interaction  of 
a  large  scale  vortical  structure  approaching  the  surface.  As  can  be  observed  from  the 
picture,  the  size  of  this  structure  is  comparable  with  the  local  width  of  the  jet.  It  is  striking 
that  similar  large  scale  structures  are  not  seen  on  the  shadowgraph  image  upstream  of  the 
interaction  with  the  free  surface.  A  plausible  explanation  is  that  the  intensity  variation  on 
the  shadowgraph  image  is  due  to  variations  of  the  second  order  spatial  derivative  of  the 
index  of  refraction  integrated  along  the  beam  path.  It  follows  that  the  shadowgraph  image 
of  the  turbulent  region  enhances  the  small  scale  structures  in  the  flow.  As  the  turbulent  jet 
reaches  the  surface  the  refraction  effects  associated  with  the  surface  defcxmation  dominate 
the  shadowgraph  image.  This  surface  deformation  is  caused  by  the  large  scales  in  the  jet 
flow.  Also  from  Figure  3.6,  a  measure  of  the  visual  growth  rate  of  the  jet  upstream  of  the 
interaction  can  be  determined  which  gives  a  value  of  1 1  degrees  for  the  half  angle. 

Figures  3.7  and  3.8  are  shadowgraph  images  of  the  flow  for  h/d=S.S  and 
Reynolds  number  of  approximately  9.5x103  and  1.6x10^  respectively.  At  the  lower  jet  exit 
velocity,  Figure  3.7  shows  surface  deformation  caused  by  the  large  scale  turbulent 
structures  in  the  jet  at  an  axial  location  of  x/d»25  on  the  left  side  of  the  photograph.  Weak 
surface  waves  are  observed  propagating  at  an  angle  of  ^2  degrees  relative  to  the 
downstream  direction.  Farther  downstream  the  vortical  fluid  in  the  large  scale  structures 
breaks  the  surface  resulting  in  complicated  surface  patterns  (x/d*33).  Again  the  dark  spots 
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associated  with  the  vortex  lines  terminating  at  the  free  smface  are  clearly  visible  in  this 
photograph.  At  higher  jet  exit  velocity,  in  Figure  3.8,  the  surface  deformations  caused  by 
the  large  scale  structure  in  the  jet  results  in  stronger  surface  waves  being  generated  at  the 
interaction.  The  waves  propagate  at  an  angle  of  54  degrees  relative  to  the  downstream 
direction.  In  this  Figure  the  wavelength  of  the  surface  waves  is  approximately  7  cm 
compared  to  values  of  the  order  of  2  cm  measured  in  the  near  field  interaction  of  Figures 
3.4  and  3.5.  The  location  where  waves  are  first  observed  x/d*25  and  the  location  where 
the  vortical  fluid  breaks  the  surface  x/d=33  are  approximately  the  same  as  for  the  low 
velocity  conditions  in  Figure  3.7. 

III.2  Surface  Curvature  Measurements 

The  surface  curvature  measurements  were  conducted  to  quantify  the  various  surface 
features  produced  by  the  interaction  of  the  jet  with  the  free  surface.  The  measurements 
were  performed  by  locating  the  photodiode  on  a  matrix  of  54  points  and  recording  the 
temporal  fluctuations  of  light  intensity  at  each  point  on  the  surface.  An  area  of  22.9  cm 
xl2.7cm  was  covered  at  2.54  cm  intervals.  A  typical  time  trace  of  the  photodiode  output  as 
recorded  by  the  Digital  Data  Acquisition  System  is  shown  in  Figure  3.9.  No  attempt  was 
made  to  calibrate  this  signal  in  terms  of  the  surface  curvature.  However  all  the  data  were 
obtained  at  the  same  conditions  of  illumination,  amplifier  gain  and  with  the  photodiode  at 
the  same  distance  to  the  surface  so  that  meaningful  comparisons  between  the  results  of 
various  tests  can  be  conducted.  The  mean  value  of  the  photodiode  output  signal,  Vl,  is 
related  to  the  local  light  intensity  of  the  collimated  beam,  Iq.  The  temporal  fluctuation  of  the 
output  is  the  result  of  temporal  fluctuations  of  the  surface  curvature  at  the  probe  location. 
Therefore  the  rms  value  of  the  photodiode  output,  Vnns>  is  ^  measure  of  the  rms  amplitude 
of  the  surface  curvature  fluctuations. 

The  rms  value  of  the  photodiode  output  fluctuation  were  mapped  over  the  region  of 
interaction  of  the  jet  with  the  free  surface  at  several  flow  conditions.  The  results  are 
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presented  in  Figures  3. 10  to  3. 16  in  the  form  of  contour  plots  of  the  rtns  value  at  each  flow 
condition.  Also  shown  in  each  plot  is  a  straight  solid  line  which  corresponds  to  the  visual 
growth  of  the  free  jet  determined  on  the  shadow  image  in  Figure  3.6.  A  common  feature  of 
all  contour  plots  is  that  along  the  centerline  the  rms  value  reaches  a  maximum  some  distance 
from  the  nozzle  and  then  decreases  farther  downstream.  This  corresponds  to  the  region  of 
maximum  surface  activity  on  the  flow  pictures.  The  lateral  extent  of  surface  activity  region 
clearly  extends  beyond  the  visual  growth  of  the  free  jet  This  is  due  to  the  propagation  of 
waves  away  from  the  interaction  region.  The  surface  area  bounded  by  a  contour  line  for 
relatively  small  values  of  the  rms  can  be  loosely  characterized  as  a  sector  of  a  circle.  The 
straight  lines  on  the  upstream  side  of  the  sector  corresponds  to  the  direction  of  propagation 
of  waves  at  the  beginning  of  the  interaction  region.  This  angle  changes  with  flow 
conditions. 

Figure  3.10  was  obtained  at  Re=3.2xl0^  and  jet  exit  depth  corresponding  to  h/d=l 
At  these  flow  conditions  the  maximum  interaction  occurs  in  the  near  field  x/d=4.  Figures 
3.1 1  to  3.13  were  obtained  at  h/d=2.5  for  jet  exit  velocities  of  100  cm/s,  150  cm/s  and  200 
cm/s  respectively.  Figures  3.14  to  3.16  were  obtained  at  h/d=3.5  for  the  same  jet  exit 
velocities.  In  all  cases  the  jet  exit  diameter  was  0.64  cm.  The  corresponding  Reynolds 
number  varied  from  6.4x10^  to  1.27x104  as  indicated.  For  a  fixed  depth,  as  the  jet  exit 
velocity  is  increased  the  location  of  the  maximum  interaction  moves  downstream  and,  also, 
the  propagation  angle  of  waves  generated  at  the  beginning  of  the  interaction  region  is 
increased.  This  behaviour  is  found  for  both  h/d=2.5  (Figures  3.1  i  to  3.13)  and  h/d=3.5 
(Figiues  3. 14  to  3. 16).  The  effect  of  increasing  depth  for  a  constant  jet  exit  velocity  is  a 
displacement  of  the  maximum  interaction  point  farther  downstream  and  to  decrease  the 
propagation  angle  of  waves  generated  at  the  beginning  of  the  interaction. 

A  better  measure  of  the  rms  value  of  the  surface  curvature  fluctuation  is  obtained  by 
normalizing  Vnns  with  the  local  mean  light  intensity,  Vl.  The  downstream  evolution  along 
the  centerline  of  the  normalized  rms  values  of  the  surface  curvature  fluctuation  Ynns/^L.  for 
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several  flow  conditions  is  shown  in  Figures  3.17,  3.18  and  3.19.  The  results  plotted  in 
Figiire  3.17  correspond  to  a  normalized  jet  exit  depth  h/d=2.5,  at  several  jet  exit  velocities 
and  for  different  jet  exit  diameters  as  indicated  in  the  figure  caption.  The  results  plotted  in 
Figure  3.18  correspond  to  h/d=3.5  and  for  several  jet  exit  velocities  and  jet  exit  diameters 
as  well.  The  results  for  Figure  3.19  are  for  h/d=1.5,  Ue=200  cm/s  and  h/d=l,  11^=50 
cm/s.  The  jet  exit  diameter  for  both  cases  was  d=0.64  cm.  Several  trends  are  readily 
identified  in  the  data.  At  a  fixed  normalized  depth,  say  for  h/d=2.5  (Figure  3.17),  the 
effect  of  increasing  jet  exit  velocity  is  to  increase  the  maximum  value  of  Vjms/VL  and  to 
displace  the  location  of  the  maximum  value  downstream.  The  effect  of  jet  exit  diameter  for 
fixed  jet  exit  velocity  and  normalized  depth  is  very  small.  At  the  larger  value  of  h/d  the 
same  trends  are  observed  for  variations  of  the  jet  exit  velocity  and  the  jet  diameter.  The 
eflect  of  the  normalized  depth  h/d  on  these  profiles  is  determined  by  comparing  the  various 
curves  in  Figures  3.17,  3.18  and  3.19.  It  is  apparent  that  increasing  h/d  results  in  a 
displacement  of  the  interaction  point  in  the  downstream  direction  and,  for  the  same  jet  exit 
velocity,  a  reduction  of  the  maximum  value  of  VnnsA^L* 

The  surface  curvature  fluctuations  were  further  characterized  with  power  spectrum 
measurements  along  the  nterline  of  the  flow  and  at  selected  locations  away  from  the 
centerline  where  waves  were  observed  in  the  pictures.  Figure  3.20  is  a  typical  plot  of  the 
frequency  power  spectrum  measured  on  the  centerline  at  x/d=24  for  Ue=150  cm/s  and 
h/d=3.5.  This  point  corresponds  to  point  1  in  Figure  3.15  in  the  region  of  maximum  rms 
fluctuation  for  this  flow  conditions.  The  vertical  axis  in  this  power  spectrum  plot  is 
normalized  with  the  square  of  the  rms  value  so  that  the  area  under  the  curve  is  unity.  There 
is  a  distinct  peak  in  the  power  spectrum  at  a  frequency  of  8  Hz.  The  corresponding 
normalized  frequency  using  jet  exit  parameters  is  fd/Ue  =  0.034  at  least  an  order  of 
magnitude  below  the  value  for  the  preferred  mode  fenr  the  free  jet  at  these  conditions. 
Similar  power  spectrum  plots  were  obtained  along  the  centerline  at  other  flow  conditions. 
From  these  plots  the  frequency  corresponding  to  the  maximum  power  spectral  density  was 
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determined.  The  results  are  presented  in  Rgure  3.21  where  this  firequency  normalized  with 
the  jet  exit  velocity  and  diameter  is  plotted  as  a  function  of  x/d  for  several  values  of  h/d. 

The  diameter  of  the  jet  was  0.64  cm  in  all  cases.  For  h/d=3.5,  solid  symbols,  the 
measured  frequency  is  approximately  constant  up  to  x/d=  40  the  farthest  location  measured 
in  this  case.  In  contrast  for  h/d=2.5  and  1.5  there  is  a  significant  drop  of  the  frequency  for 
the  larger  values  of  x/d.  There  is  also  significant  scatter  among  the  normalized  frequency 
measured  at  different  flow  condidons  suggesting  that  the  jet  exit  parameters  are  not  the 
most  appropriate  scaling  parameters  for  this  flow. 

The  power  spectrum  measured  along  the  centerline  can  also  be  compared  with  a 
similar  power  spectra  measured  away  from  the  centerline.  Figure  3.22  is  such  a  power 
spectrum  measured  at  x/d=28  and  y/d=12  for  Ue=150  cm/s  and  h/d=3.5,  the  same  flow 
conditions  as  for  the  power  spectrum  in  Figure  3.20  but  at  point  2  in  Figtire  3.15.  The 
power  spectrum  shows  three  distinct  peaks  of  frequencies  of  8  Hz,  16  Hz  and  24  Hz.  'I  ne 
tallest  peak  occurs  at  the  same  frequency  as  the  peak  in  point  1  (Figure  3. 15)  on  the 
centerline.  The  fact  that  the  frequencies  are  the  same  indicates  that  the  measured  surface 
deformations  are  the  result  of  waves  propagating  through  points  1  and  2  on  the  surface. 
This  is  of  course  consistent  with  the  results  of  the  flow  visualization  study. 

ni.3  Flow  Velocity  Measurements 
in.3.1  Free  Jet  Results 

A  systematic  study  of  the  flow  field  caused  by  the  jet  at  a  large  distance  from  the 
free  surface  was  conducted  ^o  determine  the  free  jet  behaviour  in  the  facility.  As  mentioned 
earlier  the  jet  establishes  a  recirculating  flow  in  the  free  surface  tank.  While  this  effect  is 
expected  to  be  small,  the  free  jet  data  was  measured  in  the  tank  to  confum  this  expectation 
and  to  provide  basic  data  for  comparison  with  the  free  surface  cases  discussed  below.  The 
measurements  were  conducted  at  a  jet  exit  velocity  of  200  cm/s  and  jet  exit  diameter  d=0.64 
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cm.  The  corresponding  Reynolds  number  is  Re=  1.27x10^.  The  centerline  of  the  jet  was 
positioned  24  diameters  below  the  surface  (h/d=24). 

The  mean  and  rms  value  of  the  velocity  fluctuation  was  measured  along  vertical  and 
horizontal  traverses  at  several  distances  from  the  jet  exit  plane.  The  measured  mean 
velocity  profiles  are  shown  in  Figure  3.23  for  the  vertical  profiles  (z-direction)  and  in 
Figure  3.24  for  the  horizontal  profiles  (y-direction)  respectively.  In  these  figures  the 
velocity  profiles  are  normalized  with  the  maximum  mean  value  of  the  velocity,  Um, 
measured  at  each  cross-section  on  the  centerline.  The  vertical  or  horizontal  coordinate  is 
measured  from  the  centerline  and  is  normalized  by  the  jet  exit  diameter.  In  each  figure  the 
measured  profiles  at  x/d=4,  8,  12,  16, 20, 24, 32  and  40  are  presented.  Also  shown  in 
Figures  3.23  and  3.24  are  the  definition  of  the  half  velocity  width  of  the  jet  along  the 
vertical,  L^,  and  horizontal,  Ly,  directions.  As  customary  the  half  velocity  width  is  defined 
as  the  distance  from  the  centerline  to  the  point  on  the  profile  where  the  velocity  is  one  half 
the  maximum  velocity  on  the  centerline.  It  is  apparent  that  the  velocity  profiles  measured  in 
the  horizontal  and  vertical  directions  are  very  similar  which  confirms  the  axial  symmetry  of 
the  flow  and  indicates  that  the  spatial  resolution  of  the  hot  wire  probe  in  the  horizontal  and 
vertical  directions  is  adequate  for  the  study  as  mentioned  earlier  in  Section  II.4. 

The  two  main  features  of  the  mean  velocity  field  in  the  jet  is  the  lateral  growth  of  the 
profiles  with  downstream  distance  clearly  observed  in  Figures  3.23  and  3.24  and  the 
corresponding  decay  of  the  centerline  velocity,  Um.  This  effect  is  not  immediately  obvious 
in  the  profiles  in  Figures  3.23  and  3.24  because  they  are  normalized  with  the  local 
maximum  velocity.  The  downstream  evolution  of  the  measured  maximum  value  of  the 
mean  velocity  along  the  centerline  is  shown  in  Figure  3.25.  In  this  plot  the  ratio  UeAJm  is 
plotted  as  a  function  of  x/d.  The  decay  of  Um  with  downstream  distance  implies  an  increase 
of  UeAJm  with  x/d  as  shown  in  Figure  3.25.  Also  the  linear  growth  of  Ue/Um  with  x/d 
expected  from  similarity  arguments  is  found  downstream  of  x/d=12.  A  least  squares  fit  to 
the  data  for  x/di^l2  gives 
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This  result  for  the  free  jet  is  in  good  agreement  with  the  results  of  other  investigators  (e.g. 
Wygnanski  and  Fiedler  1969,  Rajaratnam  1976,  Davis  and  Winarto  1980  and  Rajaratnam 
and  Humphries  1984).  It  indicates  that  the  finite  size  of  the  facility  does  not  result  in  a 
significant  change  in  the  decay  of  the  centerline  velocity.  Although  the  momentum  flux 
may  be  reduced  as  the  flow  develops  downstream  (Schneider  1985),  its  effect  on  Um  is 
within  the  scatter  of  the  data. 

A  quantitative  measure  of  the  lateral  growth  of  the  jet  can  be  obtained  from  the 
evolution  with  downstream  distance  of  the  half  velocity  width  of  the  velocity  profiles.  The 
results  are  shown  in  Figure  3.26  where  the  normalized  half  velocity  width  in  the  vertical 
direction,  Lz/d  (open  symbols),  and  in  the  horizontal  direction,  Ly/d  (solid  symbols),  are 
plotted  as  a  function  of  the  normalized  downstream  distance,  x/d.  As  indicated  earlier  the 
measured  widths  are  almost  the  same  in  the  horizontal  and  vertical  direction  at  all 
downstream  locations.  In  the  far  field,  the  results  show  the  expected  linear  growth  with 
downstream  distance.  A  least  squares  fit  to  the  data  for  x/d>12  gives 

^=^=0.078(^+0.97) 

This  result  is  also  in  good  agreement  with  the  measured  values  reported  in  the  literature 
(e.g.  Wygnanski  and  Fiedler  1969,  Rajaratnam  1976,  Davis  and  Winano  1980,  and 
Rajaratnam  and  Humphries  1984).  It  implies  that  the  small  recirculating  flow  in  the  free- 
surface  taitk  has  no  apparent  influence  on  the  growth  of  the  fir  et.  It  is  interesting  to 
compare  the  growth  rate  determined  f  <  m  the  half  width  of  the  velocity  profile  with  the 
visual  growth  of  the  jet  obtained  in  Figure  3.6.  The  value  of  0.078  of  the  slope  of  the  half 
velocity  width  implies  a  divergence  angle  of  4.5  degrees  which  is  approximately  one  half 
the  visual  divergence  angle  of  1 1  degrees  measured  in  Figure  3.6.  Note  that  these  angles 
are  measured  with  respect  to  the  centerline  of  the  jet 
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The  mean  centerline  velocity  and  the  half  velocity  width  evolution  indicate  that 
similarity  is  obtained  for  x/d^l2.  This  implies  that  the  mean  velocity  profiles  measured 
downstream  of  x/d=12  should  collapse  on  a  single  profile  when  the  mean  velocities  are 
normalized  by  the  maximum  velocity  in  the  profile  and  when  the  transverse  coordinates  y 
or  z  are  normalized  with  the  corresponding  half-widths  Ly  or  Lz  respectively.  The 
similarity  velocity  profiles  obtained  from  these  data  are  shown  in  Figure  3.27.  There  is 
good  collapse  of  the  data.  The  scatter  of  the  data  is  due  to  uncertainty  of  the  position 
measurement  and  calibration  of  the  hot  wire  probe. 

The  turbulent  fluctuations  of  the  velocity  were  characterized  at  each  measurement 
point  by  the  rms  value.  Also  the  power  spectrum  of  the  velocity  fluctuation  was  measured 
along  the  centerline.  Figure  3.28  is  a  plot  of  rms  value  of  the  velocity  fluctuation 
normalized  by  the  mean  velocity,  both  measured  on  the  centerline  of  the  jet,  as  a  function  of 
downstream  distance.  The  turbulence  intensity  Urms/Um  increases  to  a  value  of 
approximately  0.23  for  x/d>20.  This  value  is  in  agreement  with  the  results  of  Davis  and 
Winarto  (1980).  A  value  of  UnnsAJm=0.28  was  measured  at  large  downstream  distances 
(x/d^40)  by  Wygnanski  and  Fiedler  (1969).  Comparison  with  the  mean  velocity  evolution 
indicates  that  similarity  is  reached  farther  downstream  for  the  turbulence  intensity.  These 
results  suggest  a  value  of  x/d=32  for  the  turbulent  intensity  compared  to  x/d=12  for  the 
mean  velocity  profiles.  Davis  and  Winarto  (1980)  find  a  similar  result  while  Wygnanski 
and  Fiedler  (1969)  gives  a  value  of  x/d=40  to  reach  similarity  for  the  turbulence  intensity. 

Profiles  of  the  rms  value  of  the  velocity  fluctuation  normalized  by  the  local 
maximum  mean  velocity  on  the  centerline  as  a  function  of  transverse  coordinate  normalized 
by  the  corresponding  half  velocity  width  are  plotted  in  Bgures  3.29  to  3.31.  In  these 
figures  solid  symbols  correspond  to  profiles  measured  along  the  horizontal  direction  and 
the  open  symbols  correspond  to  the  vertical  direction.  Figure  3.29  for  x/d=4  and  8  show  a 
well  defined  minimum  of  the  turbulent  intensity  on  the  centerline.  The  minimum  is  more 
pronounced  at  x/d=4  compared  to  x/d=8.  This  is  associated  with  the  termination  of 
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potential  core  region  in  the  near  field  of  the  jet  The  two  maximum  values  on  either  side  of 
the  minimum  at  the  centerline  increase  slightly  from  x/d=4  to  x/d=8.  The  rms  velocity 
profiles  at  x/d=12  and  x/d=16  are  shown  in  Figure  3.30.  The  minimum  at  the  centerline,  if 
present,  is  within  the  scatter  of  the  measurements.  The  plot  also  shows  an  increased 
normalized  rms  value  on  the  center  region  compared  to  the  profiles  measured  closer  to  the 
jet.  At  x/d=20,  24  and  32  the  results  presented  in  Figure  3.31  show  large  scatter.  This  is 
associated  with  the  increased  uncertainty  of  the  measurements  at  these  locations  where  the 
measured  velocities  are  small.  The  basic  trends  observed  in  these  figures  are  consistent 
with  measured  rms  velocity  profiles  by  Davis  and  Winarto  (1980),  Wygnanski  and  Fiedler 
(1969)  and  others.  A  common  feature  in  these  profiles  (Figures  3.29-3.31)  is  the  good 
agreement,  within  the  uncertainty  of  the  measurements,  found  between  the  horizontal  and 
vertical  profiles.  This  is  a  expected  result  from  the  axial  symmetry  of  the  free  jet  flow. 

The  power  spectrum  of  the  velocity  fluctuations  was  measured  along  the  centerline 
of  the  jet  A  distinct  peak  is  found  where  the  power  spectral  density  is  maximum.  The 
frequency  corresponding  to  this  maximum  was  determined  as  a  function  of  downstream 
location.  This  maximum  frequency  can  be  associated  with  the  passage  frequency  of  the 
large  scale  turbulent  eddies  at  a  particular  location  in  the  flow.  The  variation  of  this 
frequency  normalized  with  jet  exit  parameters,  fd/Ue,  as  a  function  of  nramalized  distance 
downstream,  x/d,  is  shown  in  Figure  3.32.  Close  to  the  nozzle  at  x/d=4  the  value  of 
fdAJc=0.36  corresponds  to  the  preferred  mode  of  the  jet  (Gutmark  and  Ho,  1983).  Far 
downstream,  x/d=12,  the  frequency  decreases  as  x’^  as  expected  from  similarity 
arguments. 

in.3.2  Free-surface  Jet  Results 

Velocity  measurements  for  the  free-surface  jet  were  conducted  at  a  jet  exit  velocity 
of  200  cm/s  and  a  jet  exit  diameter  of  0.64  cm  which  corresponds  to  a  Reynolds  number  of 
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1.27x104,  The  jet  centerline  was  positioned  at  several  depths  below  the  surface 
corresponding  to  h/d=l,  1.5,  2.5  and  3.5. 

A  typical  time  trace  of  the  velocity  measured  at  x/d=8  for  h/d=l  is  shown  in  Figure 
3.33.  The  velocity  field  at  this  and  all  other  points  were  characterized  by  the  mean  value  of 
the  velocity,  in  this  case  157.4  cm/s,  and  the  rms  value  of  the  fluctuation,  in  this  case  27 
cm/s.  In  addition  the  characteristic  time  scales  in  the  flow  were  characterized  by  power 
spectrum  measurements  at  selected  locations  in  the  flow. 

Mean  velocity  profiles  were  measured  as  a  function  of  downstream  distance  along 
directions  perpendicular  to  the  free  surface.  Figures  3.34  to  3.37,  and  parallel  to  the  free 
surface.  Figures  3.38  to  3,41,  for  various  depths  of  the  jet.  In  all  cases  the  measured  mean 
velocity  profiles  are  normalized  by  the  maximum  velocity,  Um,  measured  at  each  cross 
section.  The  transverse  coordinates  measured  from  the  jet  centerline  are  normalized  by  the 
jet  exit  diameter  d.  Figure  3.34  shows  the  velocity  profiles  measured  at  h/d=3.5  and 
several  distances  from  the  jet  exit  plane.  The  location  of  the  free  surface  is  indicated  in  the 
figure.  As  the  flow  evolves  downstream  the  jet  reaches  the  free  surface  resulting  a  mean 
velocity  at  the  surface  different  from  zero.  For  x/d<16  the  jet  has  not  reached  the  surface 
and  the  velocity  profiles  are  very  similar  to  the  free  jet  velocity  profiles.  Downstream  of 
this  station  the  normalized  mean  velocity  at  the  point  closest  to  the  surface  increases. 
However,  at  this  depth  the  maximum  mean  velocity  is  always  found  at  the  cent^line  of  the 
jet  As  the  jet  is  positioned  closer  to  the  free  surface.  Figure  3.35  for  h/d=2.5,  the  mean 
velocity  profiles  reach  the  free  surface  somewhere  downstream  of  x/d=8.  As  in  the  case 
h/d-3.5  the  normalized  mean  velocity  at  the  point  closest  to  the  surface  increases 
downstream  of  that  point.  The  maximum  velocity  on  the  profile  is  located  on  the  centerline 
at  all  distances  documented. 

The  same  trends  are  found  for  smaller  values  of  h/d.  For  h/d=1.5  in  Figure  3.36 
the  beginning  of  the  interaction  is  found  between  x/d=4  and  x/d=8.  In  this  case  the 
maximum  velocity  at  x/d=40  no  longer  occurs  on  the  centerline  of  the  jet  but  moves  closer 
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to  the  free  surface.  For  h/d=l,  Figure  3.37,  the  beginning  of  the  interaction  occurs 
upstream  of  x/d=4.  The  maximum  mean  velocity  on  the  profiles  downstream  of  x/d=24 
move  towards  the  free  surface.  As  indicated  in  Section  II.4  the  minimum  distance  to  the 
free  surface  at  which  measurements  could  be  obtained  was  1.9  mm  which  corresponds  to 
0.3  times  the  jet  exit  diameter. 

The  mean  velocity  profiles  measured  along  the  direction  parallel  to  the  free  surface 
are  shown  in  Figures  3.38  to  3.41  for  h/d=3.5,  2.5, 1.5  and  1  respectively.  The  same 
normalization  as  for  the  vertical  profiles  was  used.  These  profiles  were  measured  on  a 
plane  containing  the  centerline  of  the  jet.  As  noted  in  the  discussion  of  Figures  3.36  and 
3.37  the  mean  velocity  on  the  centerline  is  not  the  maximum  value  measured  on  the  cross 
section  for  h/d=1.5,x/d=40  (Figure  3.40)  and  h/d=l,  x/d=32  and  40  (Figure  3.41). 
However  at  these  conditions  the  centerline  velocity  is  very  close  to  the  maximum  value  on 
the  cross  section  (UAJm^«98)  so  that  in  the  plots  this  value  on  the  centeriine  can  not  be 
differentiated  from  unity.  At  a  fixed  downstream  distance  the  measured  mean  velocity 
profiles  for  various  depths  are  very  similar  to  each  other  and  also  to  the  profiles  measured 
in  the  free  jet  (Figures  3.23  and  3.24).  A  more  detail  comparison  is  conducted  below  in 
terms  of  the  half  velocity  width  measured  on  these  profiles. 

As  for  the  free  jet,  the  maximum  mean  velocity  at  each  cross  section  decreases  with 
downstream  distance.  This  evolution  of  the  mean  centerline  velocity  is  shown  in  Figure 
3.42  for  h/d=l,  1.5, 2.5  and  3.5.  In  this  figure  the  parameter  UeAJm  is  plotted  as  a 
function  of  x/d.  Also  shown  in  this  figure  is  a  solid  line  which  corresponds  to  the  least 
squares  fit  to  the  free  jet  data  discussed  above  (See  Figure  3.25).  For  x/d^l6  the  mean 
centerline  velocity  in  the  frce-surface  jet  follows  the  same  evolution  as  the  free  jet  for  all 
values  of  h/d.  Farther  downstream  the  values  of  UeAlm  for  the  free  surface  jet  are 
generally  lower  than  the  corresponding  value  for  the  free  jet  (solid  line).  This  result  implies 
that  the  mean  centerline  velocity  decay  is  slower  in  the  free  surface  jet  compared  to  the  free 
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jet  The  effect  is  more  pronounced  for  the  smaller  values  of  h/d.  For  h/d=3.5  the  departure 
from  the  free  jet  data  can  only  be  observed  for  x/d^32. 

The  mean  velocity  profiles  showed  a  displacement  of  the  location  of  the  maximum 
mean  velocity  towards  the  free  surface.  This  effect  is  quantified  in  Figure  3.43  where  the 
distance  from  the  location  of  the  maximum  mean  velocity  to  the  free  surface  Zm  normalized 
by  the  jet  exit  diameter  is  plotted  as  a  function  of  x/d  for  the  various  cases  investigated.  As 
noted  earlier,  for  h/d=2.5  and  h/d=3.5  the  location  of  the  maximum  mean  velocity  remains 
on  the  centerline  and  therefore  Zm/d=h/d.  For  h/d=l  and  1.5  the  maximum  mean  velocity 
moves  towards  the  free  surface  starting  at  x/d=20  and  32  respectively.  Also  for  h/d=l  and 
x/d=40  the  location  of  the  maximum  mean  velocity  is  as  close  to  the  surface  as  it  was 
possible  to  measure  it 

It  was  also  noted  in  the  velocity  profiles  of  the  free-surface  jet  that  the  normalized 
mean  velocity  measured  at  the  closest  distance  to  the  free  surface  increased  with 
downstream  distance.  However,  because  of  the  normalization  with  the  maximum  velocity 
on  the  cross  section  this  result  conceals  important  features  of  the  evolution  of  the  mean 
velocity  at  the  free  surface.  A  better  characterization  of  the  mean  velocity  at  the  surface  Uj 
is  obtained  when  Us  is  normalized  by  the  jet  exit  velocity  Ue.  These  results  are  presented 
in  Figure  3.44  wnere  UsAJe  is  plotted  as  a  function  of  x/d  for  all  values  of  h/d  investigated. 
At  a  fixed  value  of  h/d  the  mean  velocity  at  the  surface  increases  to  a  maximum  value  and 
decays  further  downstream.  As  h/d  is  increased  the  maximum  value  of  Us/Ue  is  reduced 
and  the  location  of  this  maximum  value  moves  downstream. 

These  results  on  the  mean  velocity  profiles  of  the  free-surface  jet  show  that  as  the 
depth  of  the  jet  is  increased  the  interaction  region  moves  downstream,  which  in  turn  results 
in  a  reduction  of  the  maximum  mean  velocity  measured  on  the  free  surface.  This  is  in 
general  agreement  with  the  results  of  the  surface  curvature  measurements.  From  the  point 
of  view  of  the  mean  velocity  field  the  interaction  region  can  be  charactoized  by  a  reduced 
decay  rate  of  the  maximum  mean  velocity  a  downstream  displacement  of  the  location  of  this 
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tnaximum  and,  perhaps  more  precisely,  by  a  maximum  of  the  mean  velocity  at  the  free 
surface. 

The  growth  rate  of  the  free-surface  jet  was  characterized  by  the  half  velocity  width 
of  the  mean  velocity  profiles  in  the  directions  perpendicular  and  parallel  to  the  free  surface 
measured  from  the  centerline  of  the  jet.  The  definition  of  the  half  velocity  width  is  the 
same  as  the  one  used  for  the  free  jet  and  is  shown  schematically  in  Figures  3.37  and  3.41 
for  the  direcdon  perpendicular  and  parallel  to  the  free  surface  respectively.  The  measured 
half  velocity  width  normalized  by  the  jet  exit  diameter  are  plotted  as  a  function  of  x/d  in 
Figure  3.45  for  all  the  values  of  h/d  considered  in  this  study.  Also  plotted  in  this  figure  is  a 
solid  line  representing  the  least  squares  fit  to  the  free  jet  data  (Figure  3.26).  The  half 
velocity  width  in  the  direction  perpendicular  to  the  surface  is  given  by  the  open  symbols 
and  the  width  parallel  to  the  surface  is  given  by  the  solid  symbols.  The  results  show 
generally  lower  values  of  the  half  velocity  width  by  as  much  as  20%  compared  to  the  free 
jet  results.  It  should  be  noted  that  in  this  nondimensional  coordinates  the  mean  velocity 
results  discussed  above  show  that  the  interaction  with  the  surface  occurs  at  different 
locations  depending  on  h/d.  More  suitable  nondimensional  parameters  are  introduced  in 
Section  rV.l  which  help  clarify  the  evolution  of  the  half  velocity  width  in  the  free-surface 
jet  problem. 

The  measured  rms  values  of  the  velocity  fluctuation  on  the  centerline  normalized  by 
the  local  maximum  mean  velocity  are  shown  in  Figure  3.46.  Also  shown  in  this  figure  is  a 
solid  line  corresponding  to  the  free  jet  data.  Upstream  of  x/d=12  tne  results  for  the  free- 
surface  jet  are  in  good  agreement  with  the  free  jet  data  for  all  values  of  h/d.  Downstream  of 
this  point  the  free  surface  jet  values  are  lower  than  for  the  free  jet.  It  should  be  noted  that 
the  uncertainty  of  these  results  increases  with  x/d.  The  uncertainty  of  UnnsAJm  at  x/d=32 
is  estimated  as  ±0.05. 

Profiles  of  the  rms  value  of  the  velocity  fluctuation  normalized  by  the  local 
maximum  mean  velocity  for  the  various  values  of  h/d  are  presented  in  Figures  3.47  to  3.53 
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for  x/d=4,  8, 12, 16, 20, 24  and  32  respectively.  Two  plots  are  given  in  each  figure. 

Plots  labeled  (a)  are  the  profiles  in  the  direction  perpendicular  to  the  free  surface.  Plots 
labeled  (b)  are  the  profiles  in  the  direction  parallel  to  the  free  surface.  In  each  case  the 
transverse  coordinate  is  normalized  by  the  corresponding  half  velocity  width.  The  vertical 
velocity  profiles  are  limited  on  the  right  hand  side  by  the  location  of  the  free  surface  which 
is  different  for  each  value  of  h/d  and  for  each  downstream  location.  For  x/d<12  the  shape 
of  the  rms  velocity  profiles  in  the  direction  parallel  to  the  free  surface.  Figures  3.47(b)  to 
3.53(b),  are  in  general  agreement  with  the  free  jet  profiles  for  all  values  of  h/d.  The  shape 
of  profiles  in  the  direction  perpendicular  to  the  free  surface.  Figures  3.47(a)  to  3.53(a),  are 
also  in  general  agreement  with  the  free  jet  data.  However  it  should  be  noted  that  each  of  the 
vertical  profiles  is  abruptly  terminated  at  the  free  surface.  The  variations  along  the 
centerline  noted  in  regard  to  Figure  3.46  can  also  be  observed  in  these  profiles.  Also 
increased  uncertainty  with  downstream  distance  is  apparent  in  the  larger  scatter  of  the 
results  at  the  farthest  distance  from  the  jet  exit.  This  increased  scatter  is  also  found  in  the 
free  jet  data. 

The  turbulent  velocity  fluctuations  were  further  characterized  by  power  spectrum 
measurements  along  the  centerline  of  the  jet.  A  typical  power  spectrum  of  the  velocity 
fluctuation  measured  at  x/d=8  for  h/d=l  is  shown  in  Figure  3.54.  Although  the  spectrum 
is  broadband,  there  is  a  fairly  well  defined  peak  at  a  frequency  of  30  Hz.  This  frequency  is 
associated  with  the  large  scale  structures  in  the  flow.  The  frequency  for  maximum  power 
spectral  density  was  determined  as  a  function  of  downstream  distance  along  the  centerline 
for  all  values  of  h/d.  The  results  are  presented  in  Figure  3.55  in  terms  of  the 
nondimensional  frequency  fd/Ue.  Also  shown  in  this  figure  is  a  solid  line  which 
corresponds  to  the  free  jet  data.  The  results  generally  follow  the  free  jet  trend.  However 
significant  differences  exists  without  a  well  defined  trend  which  suggest  that  jet  exit 
parameters  are  not  the  most  appropriate  scaling  parameters  for  the  free  surface  jet 
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in.3.3  Wall  Jet  Results 

Experiments  were  conducted  to  detennine  the  velocity  field  in  a  wall  jet.  A  single 
flow  condition  was  investigated.  The  jet  exit  velocity  was  Ue=200  cm/s  and  the  jet  exit 
diameter  was  d=0.64  cm.  The  jet  was  positioned  at  a  distance  from  the  solid  wall 
corresponding  to  h/d=l. 

Rgure  3.56  and  3.57  show  the  profiles  of  the  local  mean  velocity  normalized  by  the 
maximum  mean  velocity  at  each  cross  section  measured  along  directions  perpendicular  and 
parallel  to  the  solid  surface  respectively.  The  horizontal  profiles  in  Figure  3.57  were 
obtained  on  a  plane  parallel  to  the  solid  surface  at  the  distance  from  it  where  the  maximum 
mean  velocity  was  measured.  In  contrast  with  the  free-surface  jet  the  vertical  profiles  show 
the  development  of  a  boundary  layer  on  the  solid  surface  so  that  the  maximum  velocity  is 
found  at  a  distance  from  surface  at  all  downstream  locations.  Note  that  the  z  coordinate  in 
this  figure  is  measured  with  respect  to  the  solid  surface.  The  profiles  parallel  to  the  wall 
have  the  same  general  shape  as  the  free  jet  and  free-surface  jet  Although  the  growth  rate  in 
this  case  is  much  larger  than  in  those  other  two  cases  as  is  immediately  evident  from  these 
profiles.  This  aspect  of  the  flow  will  be  quantified  by  the  half  velocity  width  results 
discussed  below. 

The  decay  of  the  maximum  mean  velocity ,Utn,  with  downstream  distance  for  the 
wall  jet  is  shown  in  Figure  3.58.  As  in  other  cases  the  parameter  Ue/Um  is  ploned  as  a 
function  of  x/d  in  this  figure.  Also  shown  in  the  figure  is  a  solid  line  which  corresponds  to 
the  least  squares  fit  to  the  free  jet  data  in  Figure  3.25.  The  results  for  the  solid  wall  jet  are 
in  good  agreement  with  the  results  of  Davis  and  Winarto  (1980).  It  is  apparent  that 
downstream  of  x/d=15  the  maximum  mean  velocity  for  the  wall  jet  is  larger  than  for  the 
free  jet  at  same  downstream  distance.  This  behaviour  is  similar  to  the  free-surface  jet  A 
detailed  comparison  with  the  free-surface  jet  will  be  presented  in  Section  IV.3  after  the 
proper  scaling  parameter  for  these  flows  are  introduced. 
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The  location  of  the  maximum  mean  velocity  normalized  by  the  jet  exit  diameter  as  a 
function  of  normalized  downstream  distance  x/d  are  shown  in  Rgure  3.59.  The  results 
show  a  rapid  displacement  of  the  maximum  mean  velocity  toward  the  surface  upstream  of 
x/d»15  followed  by  an  approximately  linear  growth  of  z^/d  with  x/d  downstream  of  this 
location.  This  result  is  also  in  good  agreement  with  Davis  and  Winarto(  1980).  It  should 
be  noted  that  the  minimum  distance  to  the  wall  is  found  at  the  same  downstream  location 
where  the  values  of  Um  depart  from  the  ftee  jet  results  in  Figure  3.58.  Also  the  linear 
growth  of  Zm  vith  downstream  distance  is  consistent  with  the  behaviour  of  three- 
dimensional  wall  jets  in  the  far  field  (Launder  and  Rodi,  1981). 

The  evolution  of  the  half  velocity  width  normalized  by  the  jet  exit  diameter  is  shown 
in  Figure  3.60  for  directions  perpendicular  (open  symbols)  and  parallel  (solid  symbols)  to 
the  wall.  In  this  case  the  half  velocity  width  is  measured  relative  to  the  location  of  the 
maximum  velocity.  Again  these  measurements  are  in  good  agreement  with  the  results  of 
Davis  and  Winarto  (1980).  It  is  apparent  that  the  wall  jet  grows  much  faster  in  the  direction 
parallel  to  the  wall  compared  to  the  direction  perpendicular  to  the  wall.  In  the  far  field  and 
at  sufficiently  high  Reynolds  number  these  length  scales  are  expected  to  vary  linearly  with 
downstream  distance  (Launder  and  Rodi,  1981).  The  results  show  an  approximately  linear 
behaviour  for  Lz  as  a  function  of  x  with  slope  dLz/dx  =  0.040.  The  expected  linear 
dependence  of  Ly  on  x  is  ntx  yet  found  for  x/d=40.  An  estimate  of  dLy/dx  obtained  from 
the  last  two  points  of  the  curve  gives  dLy/dx=0.30.  These  values  can  be  compared  with  the 
value  for  a  free  jet  0.078.  Thus  the  growth  rate  in  the  direction  perpendicular  to  the  wall  is 
0.5 1  times  the  value  for  the  ftee  jet  and  the  growth  rate  in  the  direction  parallel  to  the  wall  is 
3.9  times  the  free  jet  value.  These  result  are  in  good  agreement  with  the  results  of  Davis 
and  Winarto  (1980)  (dLy/dx=0.33  and  dLz/dx=0.036)  as  well  as  the  asymptotic  values 
reported  by  Launder  and  Rodi  (1981)  in  their  wall  jet  data  review  for  a  three-dimensional 
wall  jet  (dLy/dx=0.26  ±0.02,  dLj/dx  =0.039  ±  0.003). 
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The  downstream  evolution  of  the  rms  value  of  the  velocity  fluctuation  at  the 
maximum  mean  velocity  location  normalized  by  the  local  mean  velocity  is  shown  in  Figure 
3.61.  Also  plotted  in  this  figure  is  a  solid  line  corresponding  to  the  free  jet  values. 
Downstream  of  x/ds»10  the  measured  turbulent  intensity  is  lower  than  the  corresponding 
value  for  the  free  jet.  For  x/d^20  the  wall  jet  values  are  13%  lower  than  for  the  free  jet. 
Comparison  with  the  free  surface  jet  data  at  h/d=l  (Figure  3.46)  show  lower  values  of 
Urms/Um  fof  the  wall  jet  in  the  region  8<x/d<20.  Downstream  of  x/d»20  the  turbulent 
intensities  are  similar  in  both  cases.  As  noted  earlier  the  uncertainty  of  these  measurements 
at  x/d*32  is  ±0.05. 

The  profiles  of  the  rms  value  of  the  velocity  fluctuation  normalized  by  the  maximum 
mean  velocity  at  each  cross  section  arc  shown  in  Figure  3.62  for  x/d  =4  and  8,  and  in 
Figure  3.63  for  x/d=16,24  and  32.  In  each  figure  plot  (a)  represents  the  profiles  measured 
along  the  direcdon  perpendicular  to  the  wall  and  plot  (b)  represents  the  profiles  measured  in 
the  direcdon  parallel  to  the  wall.  The  profiles  in  the  direcdon  perpendicular  to  the  wall  are 
terminated  for  negadve  values  of  z/l^  at  the  wall.  For  x/d^8,  Figure  3.62,  the  measured 
values  of  turbulent  intensides  are  comparable  with  values  in  the  free  jet  or  the  free-surface 
jet.  On  the  other  hand  in  Figure  3.63  (b)  for  x/d^l6  the  rms  velocity  profile  in  the 
direcdon  parallel  to  the  wall  shows  a  disdnet  plateau  for  y/Ly  =  ±2  which  is  not  found  in 
the  free  jet  or  free-surface  jet  results.  This  difference  in  the  rms  velocity  profile  and  the 
much  larger  growth  rate  in  this  direcdon  suggest  a  fundamental  change  in  the  turbulent 
structure  of  this  flow  compared  to  the  free  jet  or  the  free-surface  jet.  This  aspect  of  the  wall 
jet  will  be  discussed  tii  Secdon  IV.3. 


CHAPTER  IV 


DISCUSSION 


Perhaps  the  most  interesting  result  of  this  investigation  is  the  complex  pattern  of 
surface  waves  and  motions  produced  by  the  underwater  jet  on  the  free  surface.  This 
surface  pattern  is  driven  by  the  underwater  flow.  The  flow  characteristics  and  scaling  of 
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this  underwater  flow  are  discussed  first  in  Section  TV.l.  This  is  followed  by  the 
discussion  of  the  free  surface  phenomena  in  Section  rv.2.  Finally  in  Section  IV.3  a 
comparison  of  the  free- surface  jet  and  the  wall  jet  is  presented  and  discussed. 

IV.l  Mean  Flow  Characteristics  and  Scaling  of  the  Free  Surface  Jet 

The  jet  flow  structure  was  altered  by  the  interaction  of  the  jet  with  the  free  surface. 
The  downstream  evolution  of  the  maximum  mean  velocity  presented  in  Figure  3.42  is 
different  from  that  of  the  free  jet  The  maximum  mean  velocity  decays  at  a  slower  rate  than 
for  the  free  jet.  The  location  of  the  maximum  mean  velocity  approaches  the  surface  in  the 
downstream  direction,  Figure  3.43,  and  eventually  the  maximum  velocity  occurs  at  the 
surface.  A  simple  model  is  proposed  based  on  dimensional  reasoning  and  similarity 
concepts  which  describes  the  scaling  in  the  far  field  of  the  free  surface  jet  A  schematic 
diagram  of  the  model  is  shown  in  Figure  4.1.  In  this  model  the  free  surface  is  assumed  to 
be  a  plane  of  symmetry  (Novikov,  1988)  for  the  flow  with  a  momentum  flux  twice  the 
nKxnentum  flux,  Jq,  of  the  corresponding  free  jet.  It  is  further  assumed  that  the  dominant 
length  scale  is  h  the  distance  from  the  jet  centeriine  to  the  free  surface.  The  jet  exit  diameter 
plays  an  indirect  role  through  its  effect  on  the  jet  iix>mentum  Jo- 
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The  siinilarity  scaling  in  the  far  field  of  the  turbulent  axisynunetric  jet  has  been 
discussed  by  several  authors  (e.g.,  Rajaratnam  1976,  Tennekes  and  Lumley  1972  and 
Townsend  1956).  If  the  jet  momentum  flux  ^constant,  the  linear  growth  of  length  scales 
with  downstream  distance  implies  that  sufficiendy  far  downstream  compared  to  the  jet  exit 
diameter  the  mean  centerline  velocity,  Um,  can  be  written  as: 


where  po  is  the  fluid  density  and  ci,  Xo  are  constants  determined  experimentally.  In  this 
investigation  the  value  of  ci=0.162  was  found  for  the  free  jet,  (Figure  3.25).  It  follows 
that  for  the  fiee  surface  jet  at  sufficiendy  large  distance  compared  to  the  jet  depth  h,  the 
maximum  velocity  Um  is  given  by 


(IV.2) 


where  the  factor  2  is  needed  to  account  for  the  momentum  of  the  image  jet  above  the 
surface.  The  constant  ci  should  be  the  same  as  for  the  free  jet  while  the  value  of  Xq 
depends  on  the  geometry  of  the  jet  and  consequendy  can  not  be  expected  to  be  the  same  as 
for  the  free  jet 

Equation  IV.2  can  be  written  as 


This  equation  is  based  on  the  assumption  that  Jo  is  a  constant  independent  of  x.  As 
discussed  by  Kotsovinos  (1976,1978)  this  fails  to  account  for  the  momentum  flux  of  the 
entrained  fluid  which  tends  to  reduce  the  momentum  flux  as  the  flow  evolves  downstream. 
Also  in  the  fiee-surface  jet  problem,  surface  waves  generated  at  the  interaction  will  carry 
momentum  away  from  the  jet  which  if  sufficiendy  large  will  result  in  a  lower  effective 
value  of  Jq.  Also  the  presence  of  surface  active  agents  may  contribute  to  a  reduced 
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momentuni  flux.  Not  withstanding  these  observations  the  similarity  argument  suggests 
that:  (i)  the  proper  velocity  scale  for  the  free-surface  axisymmetric  jet  in  the  far  field  is 
Ued/h;  (ii)  the  proper  length  scale  is  h  the  depth  of  the  jet;  and  (iii)  the  maximum  mean 
velocity  is  found  at  the  free  surface. 

Figure  4,2  is  a  plot  of  UedAJmh  as  a  function  of  x/h  for  all  the  free-surface  jet  data 
(h/d=l,  1.5, 2.5  and  3.5).  It  is  apparent  that  the  proposed  similarity  scaling  results  in  good 
collapse  of  the  data  shown  in  Figure  3.42.  The  results  for  x/h^l5  seem  to  follow  a  straight 
line.  A  least  squares  fit  to  these  data  gives  a  slope  of 0.099.  This  value  is  somewhat 
lower  (14%)  than  the  value  calculated  from  ciHl  =  0.1 15  using  ci=0.162  the  value  for  the 
free  jet  This  result  is  somewhat  striking  since  from  the  arguments  presented  above  the 
various  processes  that  can  invalidate  the  assumption  of  constant  momentum  flux  Jo  will 
result  in  the  reduction  of  momentum  flux  and  consequently  an  increase  of  the  slope  of 
UetiAJmh  vs.  x/h. 

To  further  examine  this  question  the  mean  velocity  similarity  profiles  for  various 
values  of  x/h  are  presented  in  Figures  4.3  and  4.4.  In  each  figure,  plot  (a)  presents  the 
similarity  profiles  in  the  direction  normal  to  the  free  surface  and  plot  (b)  presents  the 
profiles  parallel  to  the  free  surface.  Fot  a  normalized  distance  of  x/h»12.  Figure  4.3(a) 
shows  a  significant  reduction  of  the  mean  velocity  close  to  the  free  surface  (the  free  surface 
is  located  at  z/Lz  *1).  Gearly  this  velocity  profile  is  not  consistent  with  the  assumptions 
used  to  derive  equation  rv.3.  Note  that  x/h»12  is  where  the  maximum  mean  velocity  for 
the  free-surface  jet  begins  to  deviate  from  the  free  jet  line.  Farther  downstream  at  x/h=24 
and  32  the  mean  velocity  similarity  profiles  arc  given  in  Figure  4.4.  The  similarity  profile 
in  the  direction  perpendicular  to  the  surface.  Figure  4.4(a),  at  x/h=24  still  shows  the 
maximum  velocity  away  from  the  surface.  At  x/h=32  the  maximum  occurs  closer  to  the 
free  surface.  It  follows  that  the  scaling  given  by  equation  IV.3  can  only  be  expected  to  be 
valid  downstream  of  x/hs32.  If  in  Figure  4.2  the  last  two  points  are  used  to  estimate  the 
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slope  of  UcdAJmh  vs.  x/h  the  result  is  0.1 14.  This  value  is  in  good  agreement  with  the 
value  derived  from  the  free  jet  result  ci/V  2  =  0. 1 15  . 

The  growth  rate  of  the  mean  velocity  profile  were  characterized  by  the  half  velocity 
width  Ly  and  in  the  direction  parallel  and  perpendicular  to  the  free  surface.  This  half 
velocity  width  were  determined  with  reference  to  the  jet  centerline.  The  arguments 
presented  above  suggest  that  while  this  definition  is  adequate  for  the  direction  parallel  to  the 
free  surface,  a  more  appropriate  length  scale  in  the  direction  perpendicular  to  the  surface  is 
the  half  velocity  width  measured  from  the  free  surface  L  ^.  The  normalized  half  velocity 
widths  Ly/h  and  L  ^  /h  are  plotted  in  Figure  4.5  as  a  function  of  x/h  for  all  values  of  h/d. 
The  values  of  L  2  /h  (open  symbols)  are  along  a  line  parallel  to  the  Ly/h  data  but  displaced 
by  1.  This  is  simply  a  manifestation  of  the  fact  that  because  of  the  displacement  of  the 
centerline  with  respect  to  the  free  surface  the  jet  cross  section  is  enlarged  with  major  axis  in 
the  direction  perpendicular  to  the  free  surface.  Both  Ly  and  L'j  grow  approximately 
linearly  with  x.  The  growth  rates  dLy  /  dx  and  dL’2  /  dx  are  very  close  to  the  value  for  the 
free  jet.  However  for  x^>32  the  values  of  Ly  and  seem  to  converge  toward  each  other. 

It  is  interesting  to  note  that  the  values  of  Ly  and  L’2  show  small  yet  systematic 
deviations  from  the  linear  growth.  For  h/d=l  the  difference  (L’j  -  Ly)  is  smaller  at  x/h=12 
to  16,  increases  for  x/h»20  and  24  and  decreases  again  for  x/h>32.  A  plausible 
explanation  for  this  phenomenon  is  the  change  in  the  eccentricity  of  the  jet  cross  section 
during  its  interaction  with  the  free  surface.  A  similar  phenomenon  has  been  reported  by  Ho 
and  Gutmark  (1985)  for  the  growth  rate  of  an  elliptic  jet  Axis  switching  between  the 
major  and  minor  axes  occurred  in  the  range  of  x/d<40  due  to  the  self-induction  of  the  large 
scale  vortical  structures  in  the  flow.  Also  the  vortex  ring  experiments  of  Bemal  and  Kwon 
(1989)  show  an  oscillation  of  the  ring  eccentricity  during  the  interaction  with  a  free  surface. 

The  downstream  evolution  of  the  rms  value  of  the  velocity  fluctuation  normalized 
by  the  local  mean  velocity  at  the  interaction  point  is  shown  in  Figure  4.6.  In  this  figure 
Hed/Uimsh  is  plotted  as  a  function  of  x/h.  This  scaling  was  used  because  it  reduces  the 
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uncertainty  of  the  results.  The  straight  line  in  this  plot  represents  the  free  jet  case. 

Although  the  uncertainty  in  the  measurements  of  Unns  at  x/h=32  is  rather  large,  there  is  an 
indication  of  a  change  in  the  slope  of  the  curve  compared  to  the  free  jet  downstream  of  the 
point  x/h=24.  This  point  is  far  downstream  of  the  point,  x/h»l  1  on  Figure  4.2,  for  which 
the  maximum  mean  velocity  starts  to  deviate  from  the  free  jet  behavior. 

The  frequency  for  maximum  spectral  energy  on  the  jet  centerline  normalized  with 
the  proper  free-surface  jet  scaling  parameters,  fh^/Ued,  is  plotted  as  a  function  of  x/h  in 
Figure  4.7  for  all  values  of  h/d.  These  normalization  provides  good  collapse  of  the  data 
even  for  values  of  x/h<10.  Also  shown  in  this  figure  is  a  straight  line  with  slope  -2 
corresponding  to  the  expected  dependence  of  this  frequency  on  x  (f  «  x'2 ), 

The  downstream  evolution  of  the  normalized  mean  surface  velocity  along  the  jet 
centerline,  measured  at  a  distance  of  approximately  2  mm  below  the  surface,  is  shown  in 
Figure  4.8.  The  mean  surface  velocity  is  very  small  for  x^  ^  5.  The  mean  surface 
velocity  reaches  a  maximum  at  x/h*l  1  and  decreases  downstream  of  this  point  The  solid 
line  in  this  plot  represents  a  least  squares  curve  frt  through  the  normalized  maximum  mean 
velocity  data  presented  in  Figure  4.2.  From  Figure  4.8  it  can  be  seen  that  the  rate  of  decay 
of  the  surface  velocity  approaches  the  value  for  the  maximum  mean  velocity  in  the 
downstream  direction.  It  is  rather  interesting  to  note  that  the  rate  of  decay  of  surface 
velocity  is  much  slower  than  its  initial  rate  of  increase  in  the  axial  direction.  Also  from  the 
point  of  view  of  velocity  field  the  location  of  the  maximum  interaction  point  can  be  defined 
as  x/h»ll. 

IV.2  Surface  Waves  and  Surface  Motions 

One  of  the  main  objectives  of  this  investigation  was  to  gain  some  understanding  on 
the  nature  of  the  free  surface  motions  caused  by  the  interaction  of  turbulence  in  the 
underwater  jet  with  the  free  surface.  One  of  the  important  aspects  of  this  interaction  is  the 
generation  of  the  surface  waves.  These  waves  are  produced  by  the  large  scale  vortical 
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motion  moving  underneath  and  approaching  the  surface,  initially  deforming  it  and 
eventually  'breaking'  the  surface.  The  series  of  events  leading  to  generation  of  surface 
waves  observed  on  the  motion  pictures  of  the  flow  are  illustrated  in  Figures  3.3  to  3.8. 
Additional  evidence  can  be  found  in  the  results  presented  in  Figures  4.9, 4.10  and  Table 
4. 1 .  In  Figure  4.9  the  normalized  frequency  of  the  peak  power  spectral  density  of  the 
surface  deformation  measured  on  the  centerline  (open  symbols)  and  of  the  centerline 
velocity  fluctuations  (solid  symbols)  are  plotted  as  functions  of  x/h  for  two  typical 
conditions.  Figure  4.10  is  the  plot  of  the  normalized  rms  amplitude  of  the  surface 
curvature  plotted  as  a  function  of  x/h  for  several  flow  conditions.  Of  course  x/h  is  the 
proper  normalization  for  this  flow  as  shown  above. 

It  is  apparent  from  Figure  4.9  that  the  surface  waves  are  initially  produced  with  a 
frequency  that  corresponds  to  the  underwater  large  scale  motion  at  x/h^S.  At  this  point  the 
mean  surface  velocity  is  very  small  as  shown  by  Figure  4.8  while  the  surface  curvature 
results  in  Figure  4.10  indicate  significant  amplitude  of  the  surface  deformation.  It  follows 
that  this  initial  surface  deformation  is  produced  before  the  underwater  vonical  flow  has 
reached  the  surface.  This  type  of  interaction  has  been  investigated  by  Tryggvason  and  Wu 
(1988,1989). 

Farther  downstream  (5<  x/h  <12)  the  frequency  of  the  surface  curvature 
fluctuations  remains  approximately  constant  while  the  frequency  of  the  velocity  flucniations 
decreases.  This  result  implies  that  the  surface  deformation  is  dominated  by  waves 
generated  upstream  propagating  into  this  region.  In  this  region  the  amplitude  of  the  surface 
deformation  increases  and  then  decreases  as  shown  in  Figure  4.10.  The  normalized 
frequency  at  the  maximum  interaction  point,  f,h^/Ued  and  the  location  of  the  maximum 
interaction  point,  x„^,  determined  from  the  surface  curvature  results  are  given  in  Table 
4.1.  Also  important  in  this  region  is  the  observation  in  Figure  4.8  that  the  surface  velocity 
increases  to  a  maximum  and  then  decreases.  Thus  the  waves  propagating  along  the 
centerplane  encounter  an  increase  in  the  surface  mean  velocity  followed  by  a  reduction  of 


the  surface  mean  velocity.  This  acceleration  and  deceleration  of  the  surface  flow  implies 
straining  of  the  surface. 

In  the  absence  of  a  mean  surface  deformation 


9  Vs  « 


(IV.4) 


Therefore  in  the  region  corresponding  to  5<x/h<l  1  on  Figure  4.8,  for  which  the  surface 
flow  is  accelerating,  dUs/9x>0,  there  is  stretching  of  the  surface  in  the  axial  direction  and 
contraction  of  the  surface  in  the  lateral  direction.  Farther  downstream  when  the  surface 
flow  decelerates,  dUs/dx<0,  there  is  contraction  of  the  surface  in  the  axial  direction  and 


stretching  of  the  surface  in  the  lateral  direction. 

The  effect  of  the  surface  currents  on  the  waves  have  been  studied  by  Evans  (1955), 
Taylor  (1955),  Hughes  and  Stewart  (1961),  Longuet-Higgins  and  Stewart  (1961),  Taylor 
(1962)  and  Peregrine  (1976).  Some  of  their  results  are  summarized  by  Phillips  (1966). 
They  found  that  when  waves  propagate  over  a  surface  with  non-unifonn  currents,  the 
waves  undergo  changes  in  wavelength,  amplitude  and  direction.  The  effect  of  increased 
surface  velocity  on  the  waves  can  be  observed  on  Figure  3.5.  The  distance  between  the 
wave  crests  increases  along  the  jet  centerline.  Also  there  is  an  indication  of  a  decrease  in 
the  curvature  of  the  crests  and  wave  amplitude  along  this  direction.  It  is  proposed  that 
surface  straining  prevents  the  propagation  and  generation  of  surface  waves  on  the  centerline 
region  of  the  jet  downstream  of  x/h-*!  1. 

Farther  downstream  for  x/h  >1 1,  Figure  4.9  indicates  that  the  surface  curvature 
fluctuations  occur  at  the  same  frequency  as  the  underwater  flow  field.  In  this  region  then 
the  surface  motion  follows  the  fluctuations  of  the  underwater  turbulence.  The  amplitude  of 
the  fluctuations  monotonically  decreases  with  downstream  distance.  A  conspicuous  feature 
of  the  interaction  of  the  vortical  structures  with  the  surface  are  the  dark  spots  associated 
with  the  vortex  lines  terminating  at  the  free  surface.  Since  the  fluid  at  the  free  surface  is 
initially  irrotational,  vonicity  at  the  free  surface  must  be  the  result  of  the  interaction  of  the 
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vortical  flow  under  the  surface  with  the  free  surface.  Figure  3.3  represents  a  typical 
photograph  of  this  interaction,  in  which  underwater  vortex  lines  in  the  near  field  of  the  jet 
opens  at  the  fice  surface.  This  type  of  interaction  is  important  because  not  only  imparts 
momentum  to  the  surface  but  also  because  of  the  associated  mass  transport.  These  flow 
processes  in  the  simpler  case  of  a  vortex  ring  have  been  investigated  by  Bernal  and  Kwon 
(1989)  and  Kwon  (1989). 

The  waves  generated  by  the  underwater  flow  also  propagate  away  from  the 
centerline.  The  frequency  of  the  waves  obtained  from  the  surface  curvature  measurements 
is  approximately  the  same  as  the  frequency  obtained  at  the  point  of  maximum  rms  surface 
curvature  fluctuation  along  the  jet  centerline  as  shown  by  the  results  tabulated  in  Table  4.1. 
The  fundamental  frequency  in  the  wave  region,  fw»  and  at  the  point  of  maximum 
interaction,  f„  were  obtained  from  the  power  spectrum  of  the  surface  curvature  data.  The 
results  for  various  jet  exit  velocities  and  depths  tabulated  in  Table  4.1  show  good 
agreement  between  these  two  frequencies.  The  wavelength  of  the  waves  was  calculated 
using  the  dispersion  relation  for  the  deep  water  waves  with  surface  tension  as  well  as 
gravity  taken  into  account  (Lighthill,  1978).  The  wavelengths  obtained  were 
approximately  1-7  cm  in  the  capillary-gravity  range.  The  Strouhal  number  corresponding 
to  the  cases  for  which  the  waves  were  generated  in  the  near  field  of  the  jet  was  calculated  to 
be  St  =  fw  dAJc  =  0.38.  This  value  of  the  Strouhal  number  falls  within  the  range  of 
St=0.24-0.64  corresponding  to  the  preferred  modes  for  jets  reported  in  literature.  (Gutmark 
and  Ho,  1983).  As  the  jet  depth  is  increased  the  Strouhal  number  defined  as  St=fw  h^AJed 
is  more  appropriate.  The  values  of  this  Strouhal  number  are  given  in  Table  4.1.  They  are 
also  within  the  same  range  of  values  as  the  preferred  mode  of  the  jet. 

The  propagation  characteristics  of  the  waves  show  a  number  of  interesting  features. 
The  flow  visualization  and  frequency  results  discussed  earlier  show  that  the  waves  are 
initially  generated  by  the  large  scale  vortical  structures  in  the  flow  before  they  reached  the 
free  surface.  It  follows  that  the  direction  of  propagation  of  the  waves  is  determined  by  the 
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ratio  of  the  advection  speed  of  the  large  scale  vortical  structures  in  the  flow  to  the  speed  of 
propagation  of  surface  waves  at  the  passage  frequency  of  the  large  scale  motion  at  the 
beginning  of  the  interaction  region  (x/h»5).  Thus  for  a  Hxed  depth,  as  the  jet  momentum  is 
increased  by  increasing  the  jet  exit  velocity  for  example,  the  advection  velocity  of  the  large 
scale  structures  is  increased  which  also  results  in  an  increase  of  the  passage  frequency  of 
the  large  scale  structures  since  the  interaction  region  does  not  move  in  space. 

Consequently,  for  values  of  h  larger  than  1.74  cm  which  results  in  surface  waves  on  the 
right  side  of  the  minimum  of  the  dispersion  relation,  the  effect  of  increase  momentum  is  to 
reduce  the  wave  speed.  This  results  in  an  increase  of  the  propagation  angle  relative  to  the 
downstream  direction.  These  arguments  are  consistent  with  the  flow  visualization  results 
of  Section  ni.l  which  shows  a  wave  propagation  angle  of  approximately  60  degrees  for 
Uc  =  50  cm/s  and  h/d=l  compared  to  39  degrees  for  Ue=35  cnVs  and  h/d=l.  The  results  of 
the  surface  curvature  measurements  (Figures  3.1 1  to  3.16)  also  show  a  steepening  of  the 
waves  propagation  angle  with  strength  of  the  interaction. 

The  measurements  of  the  surface  curvature  show  an  increase  followed  by  a 
reduction  of  the  amplitude  of  the  surface  curvature  fluctuations  with  downstream  distance 
(Figure  3.17  to  3.19  and  4.10).  As  the  'strength'  of  the  interaction  is  increased,  say  by 
increasing  the  jet  exit  velocity  keeping  the  other  parameters  constant,  the  maximum 
ncnmalized  value  of  the  rms  surface  curvature  fluctuation,VnnsA^L.  increases  in  magnitude. 
The  location  of  this  maximum  x^  moves  downstream.  These  changes  are  documented  in 
Figures  4.11  and  4. 12  for  the  maximum  value  of  Vnns/VL  and  Xn,/h  respectively.  In  both 
figures  the  strength  of  the  interaction  is  charaaerized  by  the  nondimensional  parameter 
Ued/Uwh  where  Ucd^  is  the  characteristic  velocity  scale  for  the  free-surface  jet  and  Uw  is 
the  minimum  dispersion  velocity  of  capillary-gravity  waves.  The  relevance  of  Uw  arises 
because  at  the  conditions  of  the  experiments  reponed  here  both  capillary  and  gravity  effects 
are  important  in  the  dynamics  of  the  surface  waves  and  motions.  The  results  in  Figure 
4.1 1  show  the  increase  in  Vnns/VL  to  be  approximately  like  (Ucd/Uwh)^.  The  downstream 
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displacement  of  the  locaticMi  of  the  maximum  interaction  point  reaches  Xn^hslO  for 
UedAJwh  »  5.5.  It  appears  that  x„^  may  not  increase  beyond  x/h=l  1  where  the  maximum 
surface  velocity  is  found. 


IV.3  Comparisons  of  a  Free-Surface  Jet  with  a  Wall  Jet 

Rajaramam  and  Humphries  (1984)  noted  similarities  between  the  fiee-surface  jet 
results  and  the  wall  jet  results.  Also  in  the  wall  jet  literature  (see  for  example  the  reviews 
by  Launder  and  Rodi  1981, 1983)  it  is  frequently  argued  that  the  apparent  linear  growth  of 
lAJm  with  downstream  distance  is  an  indication  of  negligible  skin  friction  at  the  wall.  This 
scaling  is  confiimed  by  the  detailed  measurements  of  Davis  and  Winano  (1980)  who  found 
a  rate  of  decay  very  similar  to  a  free  jet  If  in  fact  the  skin  friction  is  negligible  it  follows 
that  the  free-surface  jet  and  wall  jet  are  indeed  very  similar  and  Rajaratnam  and  Humphries 
comparisons  are  justified.  On  the  other  hand,  there  are  reasons  to  believe  that  these  flow 
fields  are  not  that  similar  based  on  vorticity  dynamics  at  a  free  surface  compared  to  a  solid 
surface  as  noted  in  the  introduction. 

To  address  these  issues  the  scaling  arguments  presented  in  Section  IV.l  can  be 
reworked  frx  the  wall  jet.  Figure  4.13  is  a  schematic  diagram  of  the  flow  field  for  the  wall 
jet  which  incorporates  the  image  jet  to  the  other  side  of  the  wall.  The  more  important 
difference  in  comparison  with  the  free-surface  jet  diagram  in  Figure  4. 1  is  the  effect  of  skin 
friction  at  the  wall  Xw  which  results  in  a  reduction  of  the  velocity  to  zero  at  the  wall. 
Following  the  similarity  argument,  i.c.  length  scales  increase  linearly  with  downstream 
distance,  x,  equation  IV.l  can  be  written  as 


where  ci  is  determined  by  the  shape  of  the  similarity  velocity  profiles  and  the  growth  rate  of 
length  scales  in  the  flow.  In  this  case  the  momentum  flux  J  at  a  downstream  station  will 
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have  to  incorpcH^te  the  excess  momentum  flux  associated  with  conflnement  of  the  jet  by  the 
wall  which  is  characterized  by  the  image  momentum  flux  Jq,  as  well  as  the  reduction  due  to 
skin  friction  at  the  wall  tw  Thus,  the  momentum  flux  at  x,  J,  can  be  written  as 

X  oo 

J  =  2Jo-2  Jdx  Jxwdy  (IV.5) 


It  follows  from  these  arguments  that,  if  changes  in  the  constant  ci  are  neglected,  the  ratio  of 
the  slope  of  lAJm  versus  downstream  distance  x  between  the  free  jet  and  wall  jet  is  given 
by  If  “tw  can  be  neglected  as  in  the  frce-surface  jet  then  J=2Jo  and  the  slope  of 

l/Um  is  reduced  by  a  factor  of  (1/2)^^  as  already  documented  in  Section  IV.l.  The  effect 
of  tw  is  to  increase  the  slope  of  lAJm  versus  x.  As  the  flow  develops  downstream  the 
slope  will  increase  as  friction  at  the  wall  continuously  reduce  the  momentum  flux.  It 
should  be  noted  that  at  the  downstream  location  where 

X  oo 

Jo  =  2  Jdx  Jtwdy 


the  slope  of  l/Um  will  be  equal  to  the  slope  for  a  free  jet  Therefore  the  fact  that  the  slope 
of  lAJm  is  close  to  the  value  for  a  free  jet  should  not  be  construed  to  indicate  that  the  skin 
friction  is  negligible  but  to  the  contrary  it  may  be  an  indication  that  it  is  large  enough  to 
balance  the  extra  Jq  due  to  the  conflnement  of  the  jet  Although  the  local  skin  fnction  may 
be  small  compared  to  Jie  local  dynamic  pressure,  i.e  Cf«l  (Launder  and  Rodi,  1981), 
the  integrated  effect  over  the  surface  is  the  dynamically  important  parameter  which  can  be 
expected  to  be  of  the  order  of  Jo  at  some  distance  from  the  jet  exit. 

The  assumption  that  ci  is  a  constant  may  not  be  an  accurate  assumption  for  the  wall 
jet  Two  effects  must  be  taken  into  account.  First  the  boundary  layer  on  the  wall  results 
in  a  reduced  momentum  flux  associated  with  the  similarity  velocity  profiles.  The  effect  of 
this  reduced  momentum  flux  is  to  lower  the  value  of  the  constant  ci  for  the  wall  jet  and 
consequently  to  reduce  the  decay  rate  of  Uin>  Second,  the  different  growth  rates  dLy/dx 
and  dlVdx  of  the  wall  jet  compared  to  a  fm  jet  will  also  change  the  value  of  the  constant 


Cl.  The  results  in  Hgure  3.60  indicate  a  much  larger  value  of  dLy/dx  and  reduced  value  of 
dLz/dx  compared  to  the  free  jet  It  is  easily  shown  that 

For  the  wall  Jet  the  effect  of  the  growth  rate  length  scales  is  to  increase  the  value  of  the 
constant  ci  compared  to  the  free  jet  This  implies  a  faster  decay  of  Um  with  downstream 
distance. 

A  second  consequence  of  these  scaling  arguments  is  that  the  characteristic  velocity 
scale  for  the  wall  jet  is  Ued/h  and  the  characteristic  length  scale  for  the  wall  jet  is  h,  the 
same  as  for  the  free-surface  jet  Therefore  direct  comparisons  can  be  made  between  the 
free-surface  jet  results  and  the  wall  jet  results  at  h/d=l  to  test  these  ideas. 

Figure  4.14  is  a  plot  of  UedAJmh  as  a  function  of  x/h  for  the  wall  jet,  solid 
symbols,  and  the  free-surface  jet  open  symbols.  Also  shown  in  this  plot  is  a  solid  line 
which  corresponds  to  the  free  jet  data.  For  8<x/h<30  the  maximum  velocity  for  the  wall  jet 
is  larger  than  for  the  free-surface  jet.  This  is  consistent  with  the  observation  that,  although 
the  image  effect  (confinement)  is  present  in  both  cases,  the  velocity  profiles  for  the  wall  jet 
must  be  zero  at  the  wall  resulting  in  less  nramentum  flux  which  must  be  balanced  by  an 
excess  centerline  velocity.  Farther  downstream,  x/h*24,  the  slope  of  the  wall  jet  data  is 
greater  than  for  the  free-surface  jet 

Comparison  between  the  growth  rate  of  the  free-surface  jet  and  the  wall  jet  is 
presented  in  Figure  4.15.  In  this  figure  the  normalized  half  velocity  width  for  the  free- 
surface  jet  was  measured  from  the  jet  exit  centerline.  The  wall  jet  half-velocity  width  was 
measured  from  the  location  of  the  maximum  mean  velocity.  The  solid  line  in  Figure  4.15  is 
a  least  squares  frt  to  the  half  velocity  width  data  for  the  free  jet  The  growth  rate  of  the  free 
jet  and  the  free-surface  jet  are  very  similar  up  to  x/h«32.  In  contrast  the  wall  jet  growth  rate 
is  much  larger  in  the  direction  parallel  to  the  surface  (solid  symbols)  compared  to  the 
direction  perpendicular  to  the  wall.  The  wall  jet  growth  rate  in  the  transverse  direction  is 


46 


approximately  3.9  times  higher  that  of  the  free  jet  The  growth  rate  in  the  vertical  plane  is 
0.51  times  the  growth  rate  for  the  free  jet.  These  results  for  the  wall  jet  are  in  good 
agreement  with  the  results  of  Davis  and  Winarto  (1980).  The  mean  growth  rates  of  the 
free-surface  jet  have  not  changed  as  a  result  of  the  interaction  with  the  free  surface.  This 
result  is  somewhat  different  from  the  growth  rate  results  obtained  by  Rajaramam  and 
Humphries(1984).  They  found  that  the  growth  rate  in  the  vertical  direction  is 
approximately  the  same  as  for  the  wall  jet  and  the  growth  rate  in  the  transverse  direction  is 
approximately  half  the  value  for  the  wall  jet 

Significant  differences  between  the  wall  jet  results  and  the  free-surface  jet  results 
also  are  found  in  the  location  of  the  maximum  mean  velocity.  Figure  4. 16  is  a  plot  of  the 
normalized  distance  to  the  wall  or  to  the  free  surface  of  the  maximum  mean  velocity  as  a 
function  of  x/h.  In  the  wall  jet  case  the  displacement  of  the  maximum  mean  velocity  occurs 
much  closer  to  the  wall  (x/h<20)  than  for  the  free-surface  jet  In  the  far  field,  it  is  expected 
that  for  the  wall  jet  Zm  will  continue  the  linear  growth  with  x  while,  for  the  free-surface  jet 
Znt  is  expected  to  be  zero. 

The  difference  in  the  growth  rate  between  the  wall  jet  and  the  free-surface  jet  is  very 
significant  and  is  consistent  with  the  different  dynamics  of  vorticity  at  a  free  surface 
compared  to  a  solid  wall.  In  the  case  of  the  free-surface  jet  the  interaction  of  the  vortical 
underwater  flow  with  the  free  surface  results  in  vortex  lines  terminating  at  the  free  surface 
as  shown  by  the  flow  pictures.  These  vortex  lines  remain  perpendicular  to  the  main  flow 
direction.  In  contrast,  for  the  wall  jet  vortex  lines  can  not  terminate  at  the  solid  wall,  where 
the  normal  component  of  the  vorticity  vector  must  be  zero,  and  must  rearrange  to  become 
parallel  to  the  solid  surface.  This  suggests  the  existence  in  the  wall  jet  problem  of  a  strong 
secondary  flow  on  cross-sectional  planes  as  first  discussed  by  Launder  and  Rodi  (1981). 
This  secondary  flow  produces  outward  (away  from  the  center  plane)  flow  near  the  wall  and 
inward  flow  at  a  distance  from  the  wall.  The  pronounced  shoulder  of  the  measured 


turbulent  intensiQr  profiles  shown  in  Hgure  3.63  (b)  can  be  attributed  to  the  presence  of 
such  secondary  flow. 


CHAPTER  V 


CONCLUSIONS 


An  experimental  investigation  of  the  interaction  of  a  submerged  turbulent  round  jet 
with  the  firee  surface  was  conducted.  Flow  visualization,  surface  curvature  measurements 
and  hot  film  velocity  measurements  were  performed  to  study  this  problem.  The  results  of 
the  free-surface  jet  experiments  were  compared  with  the  results  of  a  wall  jet  experiment.  A 
summary  of  the  main  conclusions  of  this  investigation  is  presented  below: 

•  Surface  waves  are  produced  as  a  result  of  the  deformation  of  the  surface  caused  by 
the  vortical  structures  in  the  jet  moving  underneath  the  surface.  The  surface 
waves  propagate  at  an  angle  to  the  downstream  direction.  Propagation  in  the 
downstream  direction  is  suppressed  by  the  surface  currents  produced  by  the  jet 
flow. 

•  In  the  far  field  surface  motions  are  driven  by  the  local  underwater  turbulent 
eddies  directly  interacting  with  the  surface.  Vortex  lines  terminating  at  the  free 
surface  are  a  common  occurrence  in  this  region. 

•  The  scaling  parameters  for  the  free-surface  jet  were  determined.  They  are  the  jet 
depth,  h,  and  the  characteristic  velocity  at  the  interaction  Ued/h.  The  decay  rate  of 
maximum  mean  velocity  in  the  far  field  is  reduced  by  a  factOT  of  2^^  compared  to 
the  free  jet.  The  mean  growth  rate  of  the  free-surface  jet  is  approximately  the  same 
as  the  free  jet 
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•  The  interaction  of  the  jet  with  the  solid  wall  results  in  a  fundamentally  different  flow 
field  compared  to  the  fiee-surface  jet  case.  There  is  a  higher  growth  rate  in  the 
direction  parallel  to  the  solid  surface  and  a  lower  growth  rate  in  the  direction 
perpendicular  to  the  surface.  This  is  a  consequence  of  different  vorticity  dynamics 
at  a  solid  wall  which  results  in  the  development  of  a  strong  secondary  flow.  Also 
the  maximum  mean  velocity  decays  at  a  higher  rate  compared  to  the  free-surface  jet. 
This  is  a  consequence  of  the  reduced  momentum  flux  of  the  wall  jet  due  to  skin 
friction  at  the  wall  as  well  as  the  increased  growth  rate. 
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APPENDIX  A 


ANALYSIS  OF  SURFACE  CURVATURE  MEASUREMENTS 

A,;  Analysis 

VisuaUzadon  and  quantitative  measurements  of  the  surface  deformation  was 
obtained  by  positioning  a  photographic  plate  or  a  photodiode  above  and  close  to  the  water 
surface  illuminated  from  below  by  a  collimated  beam  of  light.  The  collimated  beam  is 
altered  by  refraction  at  the  water/air  interface.  Local  depression  and  elevation  of  the  free 
surface  has  a  lens-like  effect  on  the  light  rays,  with  depression  of  the  surface  acting  as  a 
diverging  lens  and  elevation  of  the  surface  as  a  converging  lens.  The  following  paragraphs 
describe  the  relation  between  the  light  intensity  measured  by  the  photodiode  or  recorded  on 
the  photographic  film  and  the  shape  of  the  surface. 

We  consider  a  two  dimensional  surface  deformation  as  indicated  in  Figure  A.l,  and 
use  geometrical  optics  to  analyze  this  optical  system.  The  incoming  light  beam  has 
intensity  Iq-  The  deformed  water  surface  is  described  by  the  distance  T)(x)  to  the 
undisturbed  (mean)  surface.  The  photodiode  or  photographic  plate  is  located  a  distance  s  to 
the  undisturbed  water  surface. 

We  consider  a  ray  of  light  that  intersects  the  surface  at  x.  We  wish  to  determine  the 
location,  x',  where  the  ray  will  intersect  the  film  or  photodiode  plane.  From  the  geometry 
of  the  system  we  find: 

x’«x  +  (s-Ti )  tanOd  (A.l) 

In  order  to  determine  the  deflection  angle  6d  we  apply  the  law  of  refraction: 

0d-0r-Oi  (A.2) 
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where  6i  =  a  with  tana  =  dil/dx  =  ti'(x).  The  refraction  angle  0r  is  given  by: 


na  sinOr  =  nw  sinOi  (A.3) 

where  na  and  nw  are  the  indicies  of  refraction  of  air  and  water  respectively.  Then: 

tanGi 


sinGr  =  N  sinGi  =  N 


V  1  +  tan^Gi 


(A.4) 


with  N  =  nw  /  na-  It  follows: 


tanGr  = 


N  tanGi 


V  1  -  (N2.1)tan20i 


also: 


then: 


tanGi  =  t]' 


tanGr  =■ 


Nil 


(A.5) 


(A.6) 


(A.7) 


Note  that  ti'  ^  1/(N2-1)^^.  The  equal  sign  corresponds  to  the  total  reflection  condition. 
Thus  one  can  obtain  the  relationship: 


tanGd  =  T]  (x) 


N  -  V  1*  (N2.1)tT2 


Vl-  (N2.1)tT2  +  Nti’2 


(A.8) 


and: 


X  =  X  +  {s  -  Tjfx)}  T]  (x)  H(x) 


H(x)  = 


N  -  V  1-  (N2.1)ti'2 


(A.9) 
(A.  10) 


+  Nti’2 

The  amount  of  light  reaching  the  photographic  film  or  photodiode  in  a  small  surface 
element  dx'  at  location  x'is  given  by: 

I(x')  dx' =  lo  dx  (A.  11) 

where  dx  is  the  region  in  the  incident  beam  where  all  rays  going  through  dx'  pass. 


Therefore: 


where  the  absolute  value  sign  was  added  since  both  I(x')  and  Iq  are  always  positive. 

A.2  Physical  Interpretation 

Equations  (A.l)  and  (A.  12)  provide  a  simple  description  of  this  optical  system.  A 
plot  of  the  location  of  the  light  ray  on  the  image  plane,  x',  against  its  corresponding 
location,  X,  on  the  water  surface  is  shown  in  Figure  A.2.  The  straight  line  in  this  Figure 
has  slope  unity  and  corresponds  to  an  undisturbed  water  surface.  If  there  is  a  disturbance 
on  the  surface  the  second  term  in  the  equation  (A.l)  also  contributes  and  adds  a 
deformation  to  the  straight  line  shown  in  Figure  A.2.  The  solid  curve  on  Figure  A.2 
represents  the  case  of  small  amplitude  sinusoidal  disturbances  of  the  surface.  In  this  case 
the  slope  dx'/  dx  is  always  positive  and  x  is  a  single-valued  function  of  x'. 

As  the  amplitude  of  the  disturbance  is  increased  the  second  term  in  the  equation  can 
result  in  a  negative  slope  of  x'  versus  x  in  certain  regions  as  shown  in  Figure  A.2  by  the 
dotted  curve.  In  this  case  x  is  no  longer  a  single-valued  function  of  x',  i.e.  light  rays  from 
two  or  more  points  in  the  object/incoming  beam  pass  through  the  same  point  in  the  image 
plane.  This  corresponds  to  focussing  of  the  incoming  beam  to  a  point.  Furthermore,  in 
this  case  there  are  points  for  which  dx’/  dx  =  0.  At  these  points  according  to  equation 
(A.  12)  the  light  intensity  becomes  very  large.  These  points  form  characteristic  bright  lines 
or  caustics  in  the  images  obtained  using  this  technique. 

A.3  Practical  Considerations 

In  order  to  obtain  quantitative  data  using  this  technique  the  formation  of  caustics 
and  focussing  must  be  avoided.  For  a  given  surface  deformation  these  effects  can  be 
minimized  by  reducing  the  distance  s  from  the  water  surface  to  the  image  plane.  However, 
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for  sufficiently  large  deformations  of  the  surface  the  formation  of  caustics  and  focussing 
may  not  be  avoided. 

In  the  limit  of  small  surface  deformations  ti«s  and  'n'(x)«l  we  can  linearize 
equations  (A.9)  and  (A.  10)  to  obtain: 

x' =  x +  s'n’(x)  (N-1)  +  •••  (A.13) 

and  using  equation  (A.  12)  one  obtains: 

^ - (A.14) 

•o  dx  1  +  S  n  (X)  (N-1)  + 

i.e. 

^^=1-s(N-1)ti"(x)  (A.15) 

Thus,  the  intensity  at  the  measurement  plane  contains  a  fluctuating  component 
proportional  to  the  second  derivative  of  the  surface  deformation.  This  is  the  motivation  to 
use  the  phrase  "surface  curvature  measurenxjnts".  Note  also  that  in  regions  where  T|(x)  is  a 
maximum  “n  "(x)<  0  and  I(x')  will  be  larger  than  Iq  resulting  in  a  bright  region,  and  where 
qfx)  is  a  minimum  I(x')  is  smaller  than  Iq  resulting  in  a  dark  region. 
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d  (cm)  Ue  (cm)  5i  (mm)  0  (mm)  Resi 


0.133 

0.051 

255 

100.0 

0.094 

0.036 

360 

150.0 

0.077 

0.030 

450 

0.64 

200.0 

0.066 

0.026 

513 

1.27 

50.0 

0.154 

0.059 

295 

1.27 

100.0 

0.109 

0.042 

420 

1.27 

150.0 

0.089 

0.034 

510 

2.54 

10.0 

0.487 

0.188 

188 

2.54 

25.0 

0.308 

0.119 

297 

2.54 

30.0 

0.281 

0.108 

324 

2.54 

35.0 

0.261 
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Tabic  2.2  -  Characteristics  of  the  free  shear  layer  at  the  jet  exit 
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Figure  1.1  •  Schematic  diagram  of  the  Row  geometry. 


(a) 


Figure  2.2a  -  Photograph  of  the  jet  tank. 

(a)  View  from  the  front,  (b)  View  from  the  back. 
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(b) 


Figure  2.2b  -  Photograph  of  the  jet  tank. 

(a)  View  from  the  front,  (b)  View  from  the  back. 
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Fig.  2.3  -  Schematic  diagram  of  the  free  surface  tank  with  jet  tank  inside. 
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Fig.  2.4  •  Schematic  diagram  of  the  jet  tank. 
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Fig.  2.5  •  Schematic  diagram  of  piping  system. 
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Jet  tank 


Fig.  2.6  •  Schematic  diagram  of  the  level  control  apparatus. 
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Fig.  2.7  •  Schematic  diagram  of  the  apparatus  used  for  the  wall  jet  velocity  measurements. 


Fig.  2.8  •  Schematic  diagram  of  the  apparatus  used  for  shadowgraph  flow  visualization. 


Level  Control  Box 


Fig.  2.9  -  Schematic  diagram  of  the  apparatus  used  for  LBF  flow  visualization  technique. 
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Fig.  2.1 1  -  Schematic  diagram  of  surface  curvature  measurement  technique. 


Filter 


Fig.  2.12  -Block  diagram  of  the  data  acquisition  and  storage  apparatus  for 
surface  curvature  measurements. 


Hot  Film  Probe 


Fig.  2.13  •  Block  diagram  of  the  data  acquisition  and  storage  apparatus  for 
hot  nim  velocity  measurements. 
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Fig.  2.14  -  Schematic  diagram  of  the  apparatus  used  for  hot  film  velocity  measurements. 
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Jet  tank 


Fig.  2.  IS  -  Schematic  diagram  of  the  pressure  manometer  used  for 
measurements  of  jet  exit  velocity. 


Figure  3.1-  LIF  photograph  of  submerged  jet  flow.  Ue  =  10  cm/s, 

Reynolds  number  2.5x103,  h/d=l.  Region  covered  0<x/d<10. 


Figure  3.2-  LIF  photograph  of  submerged  jet  flow.  Ug  =  35  cm/s, 

Reynolds  number  8.9x103,  h/d=l.  Region  covered  0<x/d<10. 
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Figure  3.6  -  Shadowgraph  picture  of  underwater  jet  flow.  Ue  =  1 00  cm/s, 

Reynolds  number  1.27x10^,  h/d=3.5.  Region  covered  4<x/d<24. 
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Figure  3.9  -  Photodiode  output  signal  at  x/d=20,  y/d=0. 
Ug  =  150  cm/s.  Re  =  9.5x103,  h/d=2.5. 
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Figure  3.15  -  Contour  plot  of  rms  value  of  photodiode  output. 
Uc  =  150  cm/s.  Re  =  9.5x103,  h/d=3.5. 
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Hgure  3.17  -  Downstream  evolution  of  the  rms  value  of  the  surface 
curvature  fluctuation  along  the  jet  centerline. 

Solid  symbols  d*1.27  cm.  Open  symbols 
d^.64  cm.  h/d=2.5  for  all  cases. 
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Hgure  3.18  -  Downstream  evolution  of  the  rms  value  of  the  surface 
curvature  fluctuation  along  the  jet  centerline. 

Solid  symbols  dsl.27  cm.  Open  symbols 
d=0.64  cm.  h/ds3.S  for  all  cases. 
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Figure  3.19  -  Downstream  evolution  of  the  nns  value  of  the  surface 
curvature  fluctuation  along  the  jet  centerline. 

Solid  symbol  Uc=200  cm/s,  h/d»1.5.  Open 
symbol  Uc*50  cm/s,  h/d*l. 


f  noimalized  frequency  corresponding  to 
tion  along  the  jet  centeriine.  Solid  symbols 
s  h/ds2.S.  Symbol,  x  -  Ue  =*  200  cm/s,  h/d=1.5. 


Figure  3.23  -  Velocity  profiles  in  the  direction  perpendicular  to  the  surface. 
Uc  =  200  cm/s.  Re  =  1.27x10^,  h/d=24. 
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Figure  3.25  -  Downstream  evolution  of  maximum  mean  velocity  for  the 


free  jet. _ ,  least  sqaures  fit  to  the  data  Ue  =  200  cm/s, 

h/d=24,  Re=1.27  x  10  • 


Figure  3.26  -  Downstream  evolution  of  half  velocity  width  of  the  mean 

velocity  profiles  for  the  free  jet.  Solid  circles,  width  parallel 
to  the  surface.  Open  circles,  width  perpendicular  to  the 

free  surface - -  least  squares  fit  to  the  data. 

Uc  *  200  cm/s.  Re  =  1.27  x  10^ .  h/d=24. 
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Bgure  3.29  -  Profiles  of  the  nns  value  of  the  velocity  fluctuation 
noimalized  by  the  local  maximum  mean  velocity 
on  the  jet  centerline,  Um-  Open  symbols  vertical 
direction  (z).  Closed  symbols  horizontal  direction 

(y).  Ue=200  cm/s,  Re=1.27xl0^,  h/d=24,  x/d=4,  8. 
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ngure  3.30  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity 
on  the  jet  centerline,  Um>  Open  symbols  vertical 
direction  (z).  Closed  symbols  horizontal  direction 

(y).  Uc=200  cm/s,  Re=l. 27x10^,  h/d=24,  x/d=12,  16. 
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Figure  3.31  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity 
on  the  jet  centerline,  Um.  Open  symbols  vertical 
direction  (z).  Closed  symbols  horizontal  direction 

(y).  Ue=2(X)  cm/s,  Re=1.27xl0^,  h/d=24,  x/d=20,  24,  32. 
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Figure  3.33  -  Time  trace  of  the  calibrated  hot  film  velocity  signal  along  the  jet 
centerline  at  x/d=8.  Ue  =  200  cm/s,  Umean  =  157.4  cm/s, 

Urms  =  27.25  cm/s.  Re  =  1.27x10^,  h/d=:l. 


Figure  3.34  -  Velocity  profiles  in  the  direction  perpendicular  to  the  surface. 
Uc  =  200  enVs,  Re  =  1.27xl(y»,  h/d=3.5. 
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Figure  3.35  -  Velocity  profiles  in  the  direction  perpendicular  to  the  surface. 
Uc  =  200  cm/s.  Re  =  1.27x10^,  h/d=2.5. 


Figure  3.36  -  Velocity  profiles  in  the  direction  perpendicular  to  the  surface. 
Ue  =  200  enVs.  Re  =  1.27xl(H.  h/d=1.5. 
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Figure  3.37  -  Velocity  profiles  in  the  direction  perpendicular  to  the  surface. 
Uc  =  200  cm/s.  Re  =  1.27xl(y»,  h/d=l. 


Figure  3.38  •  Velocity  profiles  in  the  direction  parallel  to  the  free  surface. 
Ue  =  200  cm/s.  Re  =  1.27x10^,  h/d=3.5. 
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Figure  3.39  -  Velocity  profiles  in  the  direction  parallel  to  the  free  surface. 
Uc  =  200  cm/s.  Re  =  1 .27xl(H,  h/d=2.5. 
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Figure  3.42  -  Downstream  evolution  of  maximum  mean  velocity  fw 
the  free-surface  jet. _ ,  least  squares  linear  fit  to  the 

free  jet  data.  Uc=200  cm/s.  Re*  1 .27x  10^. 
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Figure  3.43  -  Location  of  the  maximum  mean  velocity  as  a  function 
of  x/d.  Ue=2(X)  cm/s,  Re=  1.27x10^. 


Figure  3.44  -  Downstream  evolution  of  the  surface  velocity  normalized 
by  the  jet  exit  velocity  Uc-  Uc=200  cm/s,  Re=l. 27x10"*. 
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Figure  3.45  -  Downstream  evolution  of  the  half  velocity  width  of  the 
mean  velocity  profiles.  Solid  symbols,  width  parallel 
to  the  free  surface.  Open  symbols,  width  perpendicular 

to  the  free  surface. _ _  least  squares  fit  to  the 

data  for  h/d=24.  Ue=200  cm/s.  Re  =  1.27  x  10^. 


Rgure  3.46  -  Downstream  evolution  of  normalized  rms  value  of  the 
velocity  fluctuation  on  the  jet  centerline.  Solid  curve 

represents  the  free  data.  Uc=200  cm/s,  Re=  1.27x10 
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Figure  3.47  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  pai^lel  to  the  free  siuface. 

Uc=200  cm/s.  Rc=l. 27x10^,  x/d=4. 
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Figure  3.48  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  parallel  to  the  free  surface. 

Uc=200  cm/s.  Rc=1.27xl0^,  x/d=8. 


Figure  3.49  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  parcel  to  the  free  surface. 

Ue=200  cm/s,  Re=1.27xl0^,  x/d=12. 
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Figure  3.50  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  pai^lel  to  the  free  surface. 

116=200  cm/s,  Re=1.27xl0^,  x/d=:16. 
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Figure  3.51  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  parallel  to  the  free  surface. 

Ue=200  cm/s.  Rc=1.27xl0^,  x/d=20. 
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Figure  3.53  -  Profiles  of  the  rms  value  of  the  velocity  fluctuation 
normalized  by  the  local  maximum  mean  velocity. 
Plot  (a)  profiles  perpendicular  to  the  free  surface. 
Plot  (b)  profiles  pandlel  to  the  free  surface. 

Ue«2(X)cm/s.  Re»1.27xl0^,  x/d=32. 
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Figure  3.54  -  Power  spectrum  of  the  velocity  signal  at  x/d*8;  y/d=0 
Uc  «  200  cm/s.  Re  =  1.27x10^.  h/d=l. 


Figure  3.SS  -  Downstream  evolution  of  the  normalized  frequency  along 
the  jet  centerline.  Straight  line  corresponds  to  the  slope  of 

•2.  Ue=200  cm/s,  Rc=l. 27x10^. 


Figure  3.60  -  Downstream  evolution  of  half  velocity  width  of  the  mean 
velocity  profiles  for  the  wall  jet.  Solid  symbol,  width 
parallel  to  the  solid  surface.  Open  symbol,  width  perpen 
dicular  to  the  solid  surface. 

Uc  =  200  cm/s.  Re  =  1.27  x  10  .  h/d=l. 
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Figure  3.61  -  Downstream  evolution  of  the  rms  value  of  the  velocity 
fluctuation  for  the  wall  jet.  Solid  curve  represents  the 

free  jet  data.  Uc  =  200  cm/s.  Re  =  1.27x10^,  h/d=l. 


Figure  3.62  -  Profiles  of  the  normalized  rms  value  of  the  velocity 

fluctuation  for  the  wall  jet.  Plot  (a)  profiles  perpendicular 
to  the  wall.  Plot  (b)  profiles  parallel  to  the  wall. 

Uc  =  200  cm/s.  Re  =  1.27x10^,  h/d=l,  x/d=4,  8. 
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Figure  3.63  -  Profiles  of  the  normalized  rms  value  of  the  velocity 

fluctuation  for  the  wall  jet.  Plot  (a)  profiles  perpendicular 
to  the  wall.  Plot  (b)  profiles  parallel  to  the  wall. 

Uc  =  200  cm/s.  Re  =  1.27x10^,  h/d=l,  x/d=16,  24,  32. 


Figure  4.1  -  Far  field  scaling  of  the  free  surface  jet. 
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Figure  4.2  -  Decay  of  maximum  mean  velocity  for  the  free-surface  jet. 
Solid  line,  slope  =  0.162,  corresponds  to  the  free  jet. 
Ue=200  cm/s.  Re  =  1.27  x  10  . 


Figure  4.3  -  Mean  velocity  similarity  profiles  for  the  free-surface 
jet.  Plot  (a)  profiles  peq}endicular  to  the  free  surface. 
Plot  (b)  profiles  parallel  to  the  free  surface. 

Uc  =  200  cm/s.  Re  =  1.27x10^,  x/h=10.67,  12,  13.33. 


Figure  4.4  -  Mean  velocity  similarity  profiles  for  the  firee-surface 
jet  Plot  (a)  profiles  per^ndicular  to  the  fm  surface. 
Plot  (b)  profiles  parallel  to  the  free  surface 

Uc  =  200  cm/s.  Re  =  1.27x10^,  x/h=24,32. 


Figure  4.5  -  Downstream  evolution  of  half  velocity  width  of  the  mean 
velocity  profiles  normalized  by  the  jet  depth  h.  Solid 
symbols,  width  parallel  to  the  free  surface.  Open  symbols, 

width  perpendicular  to  the  free  surface _ _  least  squares 

linear  fit  to  the  free  jet  data.  Ue=200  cm/s.  Re  =  1 .27  x  1 0‘ 


145 


Figure  4.6  -  Downstream  evolution  of  normalized  rms  value  of  the  velocity 
fluctuation  along  the  centerline  of  the  free  surface  jet.  Solid 
line,  least  squares  fit  through  the  free  Jet  data. 

Ue  =  200  cm/s.  Re  =  1.27  x  10^ 


Figure  4.7  -  Downstream  evolution  of  normalized  frequency  along  the  jet  centerline. 

Straight  line  corresponds  to  the  slope  of  -2.  Ue  =  200  cm/s.  Re  =  1 .27x  10 
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Figure  4.9  -  Downstream  evolution  of  normalized  frequencies  along  the  jet 
centerline.  Solid  symbols,  frequency  corresponding  to  velocity 
data.  Open  symbols,  frequency  corresponding  to  surface  curvature 
data.  Ue  =  200  cm/s.  Re  =  1.27  x  10 
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Figure  4, 10  -  Downstream  evolution  of  the  rms  value  of  the  surface 
curvanu-e  fluctuation  along  the  jet  centerline.  Solid 
symbols  d=  1.27  cm.  Open  symbols  d=0.64  cm. 
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Figure  4. 1 1  -  Maximum  rms  surface  curvature  fluctuation  as  a  function  of 
the  interaction  'strength'.  Solid  symbols  d=1.27  cm.  Open 
symbols  d=0.64  cm.  _ ,  slope=2. 
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Bgure  4.12  -  Location  of  the  maximum  ims  surface  curvature  fluctuation 

■ 

as  a  function  of  the  interaction  'strength'.  Solid  symbols 
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d=1.27  cm.  Open  symbols  d=0.64  cm. 

Figure  4. 1 3  -  Schematic  diagram  of  the  flow  field  for  the  wall  jet. 
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Figure  4. 14  -  Decay  of  maximum  mean  velocity  for  the  wall  jet  and  the 

firee-surface  jet.  Solid  symbol,  wall  jet.  Open  symbols,  frec- 
surface  jet  Solid  line  corresponds  to  the  free  jet 

Slopc=0.162.  Uc  *  200  cm/s.  Re  =  1.27x10^,  h/d=l. 


Figure  4.15  -  Downstream  evolution  of  half  velocity  width  of  the  mean 
velocity  profiles  for  the  wall  jet  and  the  free-surface  jet 
Solid  symbols  width  parallel  to  the  surface.  Open  symbols 
width  perpendicular  to  the  surface. 

Ue  =  200  cm/s.  Re  =  1.27x10^,  h/d=l. 


Figure  A.1  -  Geometry  of  the  deformed  water  surface. 
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■  Figure  A.2  -  Plot  of  the  location  at  which  the  ray  of  light  intersects 

1  the  photodiode  plane,  x',  versus  the  location, x,  at  which 

the  incoming  ray  intersects  the  undisturbed  surface. 
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